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Abstract 


India has experienced a large number of earthquakes in the past. The number of 
strong motion records from these earthquakes is very small, and is not sufficient to be 
used for the development of attenuation relationships and to evaluate seismic risk for a 
future earthquake. On the contrary, intensity maps are available for a substantial number 
of earthquakes and it is prudent to use these for the development of attenuation 
relationships 

In this study, 59 earthquakes from different parts of the country have been 
considered. The isoscismals are collected from different sources. All the isoseismals 
along with their intensity distance data are presented in a tabular fonn. The Indian 
subcontinent is divided into six attenuation provinces, viz.. Northwest, North, Northeast, 
Peninsular, Foredeep and Subduction. Different attenuation models are tried and best 
model has been found. For the first time, Indian subcontinent models relating magnitude 
and maximum intensity have been developed. 

Intensity Attenuation relationships have been developed for each region. A single 
equation is also developed for entire Indian subcontinent. Epicentral as well as 
hypocentral distance has been used and relationships are developed along the major and 
minor axes of the fault in addition to that for average attenuation (considering average 
radius of intensity contour as distance parameter). The derived relationships are 
compared with that of India and other countries. 

The approximate decay pattern of intensity for different provinces and distance up 
to which the intensity decay curve follows a particular trend has been found. It has been 
shown, whether it is prudent to use the relationships developed in other countries for 
seismic risk analysis in India. 
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Chapter 1 
Introduction 


1.1 General 

Seismic risk analysis for important projects like nuclear power plants, important 
bridges, and high dams, requires attenuation relationships which relate seismological data 
with parameters of engineering interest, such as peak ground acceleration, velocity, or 
response spectral ordinates. This implies that instmmental records of ground motions 
covering reasonable time spans be made available for the region of interest. However, for 
India enough strong motion database is not available and hence earthquake intensity maps 
of past earthquakes have to be used. The attenuation relations of intensity with distance 
can be obtained (e.g., Howell and Schultz, 1975), and combining those with ground 
motion to intensity correlation from regions with abundant data, one can estimate the 
seismic ground motion parameters for the region of interest {e.g., Rao et. ciL, 1982). 

Empirical intensity data from past earthquakes have often been used to predict 
ground motion for a given future earthquake with similar size and location [Malone and 
Bor, 1979] Intensity data have also been widely used to detemiine the location of 
earthquake, which occurred prior to modem seismological instrumentation (e.g., Malone 
and Bor, 1979). The intensity distribution can also be used to estimate focal depths, near 
field energy absorption coefficient, and other related parameters to assess the size of an 
event (e.g., Ambraseys, 1985). Very little research has been conducted in India. For 
seismic risk analysis of major projects like nuclear plants, important bridges and high 
dams, attenuation relations developed in other places like North America and California 
are still used. 

Many factors affect the attenuation of seismic intensity with distance. These 
include magnitude of earthquake, depth of focus, epicentral distance, rupture model of 
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source faull, topography and condition of the medium in which seismic waves propagate 
[Maugeri and Motta, 1993]. Due to these factors the earthquakes have regional 
characteristics and the attenuation laws of seismic intensity are different for different 
regions. Unless two regions have similar conditions of geology, tectonics, topography, 
earthquake source etc., one may not be able to apply the attenuation law obtained in one 
region to another region [Maugeri and Motta, 1993]. 

1.2 Methods of Describing Seismic Intensity Law 

The attenuation law of seismic intensity depends on many factors, amongst which 

the magnitude of earthquake M (or epicentral intensity /„) and hypocentral distance A (or 

the epicentral distances R) are most important. Accordingly, many expressions describing 
attenuation laws of seismic intensity have been developed. The following three are most 
frequently used: 

1 . If the strength of the earthquake is represented by the magnitude M, and assuming that 
there is a linear relation between the seismic intensity / and the magnitude M, and 
linear and logarithmic relation with the of distance R, the attenuation law of seismic 
intensity can be written in the following fonn, 

I =a + hM +c R + (l \n{R + D) (1.1) 

where a, h, and c are regression constants and D is a constant concerning the size of 
the earthquake source. In general, the term “ci?” can be neglected and the 
formulation of the above can be written as, 

I = a + hM +cl Iii{R + D) (1.2) 

2. Strength of an earthquake can be represented by the “imaginary” intensity /„ at the 
earthquake epicenter instead of the magnitude M. In this case, the attenuation law of 
seismic intensity can be written as follows. 


I = a + b p + c ln{R + D) 


(1.3) 
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where /„ is the imaginary intensity (higher than the maximum assigned intensity) at 

the epicenter of the earthquake, and Z) is a suitably assumed constant concerning size 
of the earthquake source. 

3. Equations 1.2 and 1 3 show that the decay of seismic intensity is independent of 
strength/,,. However, it has been observed that the rate of attenuation of intensity 


with distance is faster for small earthquakes than that for large earthquakes. 
Considering this, a different relationship for seismic intensity attenuation is proposed 
[Maugeri and Motta, 1993]. 


/ ' 

In 




Vo {Ri-^)_ 


where /„ is the intensity of the epicenter area, is the equivalent radius of the 
highest mapped isoseismal line, R is the equivalent radius corresponding to intensity /. 
and (/■/■. f//„are two constant obtained as follows. 


(// = 


R-k ■ Rk-i 


K>2 


and 


(1.5) 


¥ 


Rj - R, 


( 1 . 6 ) 


where is the equivalent radius of the isoseismal line (the boundary line of the 
area with intensity I ^ ). 

In all the three methods, R can be the epicentral distance or the hypocentral 

distance zl = v/ /? ' + /; " , where h is the depth of the earthquake focus. 

1.3 Scope of the Study 

In this study, earthquake data from past 59 Indian earthquakes have collected from 
many sources and different attenuation relationships are developed for Indian 
subcontinent. Keeping in mind geological conditions, the entire Indian subcontinent is 
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divided into six attenuation provinces, namely Northeast, Northern, Northwest, 
Peninsular, Foredeep, and Subduction province. 

Adopting the deterministic method and one step iterative regression analysis, 
attenuation equations are developed for each region. Different equations are developed 
involving magnitude, epicentral intensity, maximum intensity and the combination of 
them to airive at the best- fit equation. Based on these, attenuation characteristic of 
different regions have been compared. Important results are obtained like, how much 
would be the maximum damaged area if a particular strength of earthquake strikes the 
concerned place, at what distances the intensity drops by one unit at a particular place, 
whether intensity attenuation in India is similar to that of United States or not and which 
type of attenuation equation is best suited for a particular province. 


• •• 



Chapter 2 

Review of Existing Attenuation Relationships 


2.1 Introduction 

This chapter is aimed at bringing out the research carried out in the past on 
intensity attenuation. The attenuation relationships found in the literature have been 
reviewed. There are several attenuation relationships, showing intensity attenuation with 
the epicentrai distance and magnitude. Many researchers have tried to correlate the 
epicentral intensity with earthquake parameters such as focal depth, magnitude [Kaila and 
Sarkar, 1978]. Some researchers have also tried to consider the effect of local geological 
conditions {e.g., Lee and Trifunac, 1985). The basic equation is found to be in the fomi of 
either I-l^=a + bR + clogjg{R)orl = a + bM + cR + dlogig{R), where I is the 

intensity at epicentral distance R, 1^ is the epicentral intensity and a, b and c are the 

regression constants. The intensity attenuation relationships developed for India and 
United States in the past are presented in Tables 2.1 and 2.2. 

2.2 Factors Affecting the Intensity Attenuation Relationship 
2.2.1 Focal Depth 

The extent of the affected ai'eas and shapes of isoseismals depend on the focal 
depth [Rao, et. ciL, 1982]. For any earthquake the directional variation of the intensity is 
due to the variation in the sub-surface geology, whereas change in decay rate from one 
earthquake to another is likely to be due to the depth of focus. The greater the focal depth, 
slower is the fall of intensity with distance [Kaila and Sarkar, 1978]. For the same size or 
energy of earthquakes, the epicentral intensity {I depends on the focal depth of the 

earthquake and high magnitude (high epicentral intensity) is generally correlated with 
small focal depth [Burton et. ctl, 1985]. 
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2.2.2 Magnitude 

Magnitude is a measure of the size of earthquake. It is measured from amplitude 
of the earthquake waves and it represents the energy released at the source. A smaller 
magnitude earthquake can have same epicentral intensity as that of larger one, but then it 
will have smaller felt area depending upon the source dimension and the focal depth. 

2.2.3 Distance 

The variable that is generally used to characterize the propagation parameter is 
distance. The intensity or damage potential of an earthquake attenuates with distance 
from the macroseismic hypocenter or focus. The earthquake rupture can be extended over 
tens to hundred of kilometers, and a number of distance measurements have come into 
use. The various distance measures used in the literature are (Figure 2.4): (a) epicentral 
distance, R, (b) hypocentral distance. A, (c) distance to zone of energy release, R-, (cl) 
closest distance to rupture, Rn,p, and (e) closest distance to surface projection of rupture, 
Rjb. In literature, the intensity attenuation relationships are expressed in terms of either 
epicentral distance or hypocentral distance. In deriving the intensity attenuation 
relationships, it is assumed that the seismic energy radiates from a point source and hence 
hypocentral distance should be used instead of the epicentral distance. The epicentral 
distance can be used as an approximation for hypocentral distance for shallow focus 
earthquakes whose fault length is not large. 

Earthquakes whose fault length is large will have large difference in attenuation 
along the transverse and longitudinal directions since energy is radiated fi'om a line 
source instead of a point source. The difference in the value of intensity drop can be as 
high as one unit [Chandra, 1982(a)]. In the present work, epicentral as well as hypocentral 
distance is used and attenuation is expressed in both longitudinal and transverse directions 


of the fault. 
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2.3 Medium of Propagation 

The amplitude of ground motion depends on both the properties and configuration 
of the near surface strata through which seismic waves propagate. The properties having 
maximum effect on the ground motion are impedence and absorption. Impedence is the 
resistance to the particle motion and is defined as the product of density, seismic wave 
velocity, and the cosine of the angle of incidence [Reiter, 1991], Particle velocity is 
inversely proportional to the square root of impedence. As seismic waves propagate 
through a medium of increasing impedence, the increase in resistance to motion decreases 
the particle velocity and amplitude of seismic wave to preserve the energy. Keeping other 
factors constant, effect of an earthquake will be more on a low-density medium than on a 
high-density one [Reiter, 1991], 

Absorption is also tenned as damping or inelastic attenuation, which tends to be 
substantially greater for soft soils than for hard rocks [Chandra, 1982(a)]. At higher 
frequencies the impact of absorption is found to be very severe while at low frequencies it 
is less. These effects appear to balance out at frequencies dominating the peak ground 
accelerations and most empirical studies show that there is very little difference in peak 
ground acceleration recorded on rock or soil. If ground motion is dominated by high 
frequencies, such as in small earthquakes, the lower absorption of hard rock can allow 
higher peaks than other types of site [Reiter, 1991]. 

2.3.1 Local Geological Conditions 

The distortion of isoseismal shape from concentric circles is mainly due to the 
local geological conditions. The character of the upper few hundred meters of strata 
beneath any site may, under certain conditions, have a more profound effect on the level 
of earthquake shaking than the effect of earthquake magnitude. The different intensity 
scales make little or no effort to take the local geological conditions into account [Reiter, 
1991]. An examination of Modified Mercalli intensity scale shows that, considerable 
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damage to building could define intensity VII, VIII and IX, depending upon whether the 
buildings were “poorly built or badly designed”, and “ordinary substantial buildings”, or 
“structures built especially to withstand earthquakes”. No mention is made, however, of 
whether these structures are on hard rock, deep sediments, or loosely consolidated sands. 

Many investigators {e.g., Malone and Bor, 1979) have estimated generic 
corrections in expected intensity as a function of local site conditions. Evemden and 
Thomson (1985) used specific corrections to an average intensity to correlate with surface 
geology (taken from Reiter, 1991), as it appears on the geological map of California 
(Table 2.3). From the studies of Tilford et. al. (1985), it was concluded that sites located 
on the soil foundations appear to be subjected to an intensity of ground shaking about one 
unit higher on the MM scale than the sites located on rock foundations (Table 2.3). 

2,3.2 Source Dimension 

In deriving the intensity attenuation relationship it is assumed that the energy is 
radiated from a point source through a space of simple geometry (e.g., a unifonn 
hemisphere or layers) [Howell and Schultz, 1975]. Such type of source under uniform and 
isotropic soil conditions gives rise to an isoseismal pattern in the form of concentric 
circles. It has been often observed that isoseismals are elliptical in nature with their 
longitudinal axis along the length of the fault. To accommodate this fact different 
attenuation rates can be assumed along the major and minor axes of the elliptical 
isoseismals (e.g., Tilford et. al., 1985). 

An earthquake with a small source dimension may have a large observed or 
mapped epicentral intensity, although the total felt area might be rather small. On the 
other hand, an earthquake with a large source dimension may have the same mapped 
epicentral intensity but a much larger felt area. The earthquake in later case is of larger 
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magnitude. The calculated epicentral intensity (/^) values represent the strength of an 

equivalent point source. 

2.4 Ste ps in Developing Intensity Attenuation Relationships 

2.4.1 Collection of Intensity Maps 

In deriving the intensity attenuation relationship the isoseismal maps of the past 
earthquake are required. In India most of the isoseismal maps are available with the 
Geological Survey of India in Memoirs of GSI. The other sources include Indian 
Meteorological Department, National Geophysical Research Institute (Hyderabad) and 
many Indian journals. 

2.4.2 Preparation of Database 

The pattern of isoseismals would be concentric circles, with highest intensity near 
the epicenter, if the earthquake energy is released from a point source and propagates 
through a uniform isotropic soil. Since these conditions are not satisfied in reality, the real 
patterns of the isoseismals are irregular. The preparation of database is entirely dependent 
on the pattern of isoseismal maps, i.e., either elliptical or irregular. To prepare the 
database one can convert the mapped isoseismal areas into equivalent circles and obtain 
their radius for the corresponding intensity. Many researchers have adopted this approach 
because of its simplicity. 

Sometimes the trend of isoseismals is elliptical, with longitudinal axis along the 
length of the fault. In this case many researchers {e.g., Tilford and Chandra, 1985) have 
assumed different attenuation along major and minor axes, and data points were obtained 
in both the directions along and normal to the fault rupture. In case of irregular shapes, 
data points can be read along any number of desired directional lines and then averaged to 
get the average attenuation {e.g., Chandra, 1979). 
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2.4.3 Norm^ali 2 ation of Intensity Scale 

Over the years, various researchers have developed different intensity scales. 
Hence, for comparison and uniformity, the different scales have to be converted to one 
scale. In this study, the various intensities are converted to one scale using the conversion 
available in the literature [Reiter, 1991]. 

2.4.4 Division into Different Attenuation Provinces 

There are many factors affecting the intensity-attenuation relations, such as 

strength of earthquake, focal depth, and epicentral distance, rupture model of source fault, 
topography and medium of propagation. However, earthquakes have regional character. 
Unless two regions have similar conditions of geology, tectonics, topography, earthquake 
source etc., the attenuation law obtained in one region cannot be applied to another region 
[Maugeri and Motta, 1993]. Hence, the concerned area has to be divided into the different 
zones so that a particular equation is adequate to describe the attenuation relationship. In 
present study the Indian subcontinent is divided into six attenuation provinces, viz.. 
Northeast, Northern, Northwest, Peninsular, Foredeep, and Subduction province. 

2.4.5 Selection of Suitable Attenuation Relationship 

Many attenuation relationships are available in the literature as described in the 

Section 2.6. There are basically three methods of describing the intensity attenuation as 
described in the Section 1.2. Suitable attenuation relationships should be selected 
depending on the region and the adequacy of the relationship to the region. 

2.4.6 Set Initial Epicentral Intensity (Jo) Values 

Intensity attenuation relationship is more sensitive to the epicentral intensity rather 

than the maximum intensity. Many studies have shown that the assumed initial high value 
of epicentral intensity (/„) will result in a relationship with high rate of intensity 

attenuation relation and vice versa {e.g.. Burton et. al, 1985, Lee and Trifunac, 1985). 
Hence the initial values of epicentral intensity (/J should be selected appropriately. 
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There are two methods to airive at the appropriate epicentral intensity (/^) values. These 
are, 

1. Graphical Method: Following the parallelism of the intensity lines for different 
earthquakes, a smooth curve is drawn, using this curve as the reference, for a 
particular earthquake, /„ at each isoseismal intensity point was estimated and then 

averaged. Most of the authors use this method because this give quite good 
approximates. 

Figure 2.1 shows intensity-distance curves for average distance for the three 
earthquakes. Following the general parallelism of the intensity lines for different 
earthquakes, a smooth curve has been shown by dashed line. The intensity at zero 
distance for the dashed curve is 8.7. Using the dashed curve as a reference, for a 
particular earthquake, /„ at each isoseismal, intensity point was estimated and then 

averaged. For example, the five points for the earthquake event 1, in order of 
increasing distance occur at 0.50, 0.58, 0.23, 0.43 and 0.43 intensity units above the 
intensity 8.7 (dashed) curve. Therefore, from these three points, 1^ for earthquake 
event 1 will be inferred as 9.20, 9.28, 8.93, 9.13 and 9.13, respectively, and the 
average /„ is be estimated as 9.13. Similarly /„ for the other events are calculated. 

2. Empirical Method: If the magnitude of the earthquake is accurately known (from 
reliable source), one can use the epicentral intensity-magnitude relationship to get the 
epicentral intensity. Due to non-availability of instrumental records of past 
earthquakes, it is difficult to rely on the magnitude. Hence this method is scarcely 
used. 

2.4.7 Regression Analysis 

After setting the initial value of the epicentral intensity, iterative regression 
analysis has to be performed between drop of intensity (I ^ - 1) values and the 
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corresponding distance {A or R). Using the regression expression, the epicentral intensity 
for each data point of an earthquake is found and averaged to get the new epicentral 
intensity. This procedure is repeated until no significant improvement is observed for the 
drop of intensity (/^ - /) values. 

2.4.8 Method of Regression Analysis 

In literature many methods of regression analysis are available such as method of 
least square, two step regression analysis, and method of maximum likelihood. One-step 
method of least square is generally used {e.g., Gupta et. al, 1999, Chandra, 1979, Yarar 
el. al. 1984). 

2.4.8. 1 Method of Least Square 

n 

This procedure minimizes the sum of the squared error^w, (y, -y, )' , where y, 

1=1 

is the predicted value ofy, , and w, is the weight assigned to each observation y, . 

2. 4. 8. 2 Validity of Results 

2. 4. 8. 2. 1 Analysis of Residuals 

Different methods for regression analyses are derived on the assumption that the 
random eiTor term (£■/) is normally distributed. This error term is represented by the 
residuals obtained after regression. The nonnality assumption of the residual is checked 
by yf (chi-square) test. The adequacy of the model with respect to the fitted variables is 
then checked by plotting the residuals versus the independent variables {f, R). If these 

plots do not show any apparent trend, the model can be assumed adequate [Roshan, 

2000 ]. 

2. 4. 8. 2. 2 Multiple Correlation Coefficient 

The adequacy of the fit can be evaluated by the multiple correlation coefficients. 

The expression for coixelation coefficient is {e.g., Timiovska, 1992), 
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^ obs )^(Kst ^e^t) 


( 2 . 1 ) 


where is the observeij intensity, I is the mean of the observed intensity, 1^,^, is the 
estimated intensity level, and is the mean of the estimated intensity level. The value 
of the coiTelation coefficient should be in between 0 to 1 [Timiovska, 1992], 


2 4.8 2.3 Uncertainty in Prediction 

In addition to median intensity estimates, the uncertainty (standard deviation) is 
also important for seismic risk analyses. The uncertainty in the prediction is represented 
with the help of standard deviation (cr) of the predicted values. The predicted dependence 
of standard deviation is different for different models, since it depends on the data set. 

2.5 Difficulties in Deriving the Attenuation Relationships 

2.5.1 Use of Different Intensity Scale 

Over the years several intensity scales like Oldham, RF, MM and MSK were used. 
To develop intensity attenuation relationships, one has to convert the available isoseists of 
different intensity scales to one scale. Some difficulty arises in assigning the maximum 
intensity. For example, intensity VIII on RF scale is equivalent to VII on MM scale, but 
for an earthquake of maximum intensity VIII on RF scale, does.not mean that it will have 
maximum intensity of VII on MM scale. In such case the maximum intensity reaches VIII 
on MM scale on a separate isoseist to be drawn inside the innermost isoseismal. Hence 
careful examination is needed while converting intensities from one intensity scale to 
another. 

2.5.2 Maximum Intensity 

The empirical equations found in literature, show that intensity-attenuation 
relationship is quite sensitive to the maximum intensity. The maximum intensity is 
distributed on a relatively small area, and is more sensitive to population and building 



density than estimate of smaller intensities at greater distance. Since the intensity scale is 
not finely divided and does not take all factors in account, the maximum intensity may be 
mcoiTectly estimated [Anderson, 1978]. 

2.5.3 Use of Intensity as a Whole Number 

The available intensity scales designate different levels of shaking with whole 

numbers. Since the estimates of intensity near the epicenter of an earthquake are sensitive 
to population and building density, there is a high probability that the epicentral intensity 
can be inappropriately estimated [Anderson, 1978], Sometimes ambiguity in assigning the 
intensity level is dealt by assigning a dual intensity level, such as MMI VII- VIII (7.5), or 
VII. V (7.5). However, this type of practice is very rare and whether its adequacy is 
debatable [Reiter, 1991]. 

2.5.4 Source Dimension 

In deriving the intensity attenuation relationship it is assumed that the energy is 
radiated from a point source through a space of simple geometry (e.g., a unifonn 
hemisphere or layers). The equations do not give good results where distance is not large 
compared to the source dimension [Howell and Schultz, 1975]. An alternative approach 
postulates that the intensity-distance relation is the sum of two attenuation rates, one at 
short distance, the other at large distances [Kaila and Sarkar, 1978]. Some of the 
researchers have assumed that the different attenuation rate exists along the major and 
minor axes of the isoseismals [Tilford and Chandra, 1985]. 

2.5.5 Initial Epicentral Intensity (/<,) 

The level of the attenuation curve is quite sensitive to the initial epicentral 
intensity set selected for the analysis. The empirical relationship shows a high rate of 
attenuation if the initial epicentral intensity (/J set contains high values for different 
earthquakes, and low rate of attenuation if the initial epicentral intensity set contains low 
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values [Chandra, 1979]. Hence to avoid such ambiguity, a graphical method, as described 
in the Section 2.2.7 has been adopted by many researches to arrive at the initial set of . 

2.6 Literature Review of .Attenuation Relationships 

2.6.1 Ergin (1969) 

The observed intensity-distance relation for a number of earthquakes that occurred 
in Turkey and elsewhere was studied and found that the curve, given by the 

"yj ” I j 

fomiula/ - 7,^ = /; log n, -i- absorption term , best fits for either n equal to 3 or 5, 

b 

where h is the focal depth. ' 

2.6.2 Rutlage (1972) 

The macroseismic areas for all earthquakes of MM Intensity IV or greater which 
occuired in the United States (east of the longitude 106° W) from 1928 to 1969 have been 
examined. A regression equation between the epicentral distance, R, and the epicentral 
intensity, /„ , was fit. The equation used was of the form, 

= + ( 2 . 2 ) 

2.6.3 Gupta and Nuttli (1976) 

A deterministic approach was adopted to derive the intensity attenuation 
relationship for central United States by considering eight earthquakes, and adopting the 
method of least square for the regression analysis. The equation derived was, 

1-1^=3.7-0.0011 R-2.7 log, o{r) for7?>20km. (2.3) 

This equation was also used in estimating the topographical error of an earthquake whose 
maximum intensity is not known reliably. 

2.6.4 Anderson (1978) 

A probabilistic approach was adopted to develop the attenuation of shaking 
intensity with distance for the United States earthquakes. The best fitting parameters were 
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found for a Gaussian distribution function to fit the observed distributions of log,fj{R) for 

each of combination of epicentral intensity and the intensity bounded by the counters. 
2.6.5 Kaila and Sarkar (1978) 

A comprehensive atlas of isoseismal maps of 26 major earthquakes in India was 
prepared. These isoseismals were also superimposed on the geological map of India. 
Intensity-distance relations for these earthquakes for different directions were studied. 
Assuming that the logarithm of intensity varies linearly with the epicentral distance a 
single attenuation relationship for entire India was derived. The equation was, 


1 

— = exp 



+ 0 0011 J 


for 15</?<240 km, 


(2.4a) 


where /„ is the epicentral intensity in MM, I is intensity at an epicentral distance A, and 

h is the focal depth in km. The derived empirical relation connecting maximum intensity, 
/„ , magnitude, M, and focal depth, h, of an earthquake was as follows, 

I„=}.5M-4.1log {h)+4.5 for 8 < A < 70 km. (2.4b) 

A linear relationship was proposed between magnitude and inverse of the square of focal 
depth was, 

-C = -1.1 706 M -0.005. ' (2.4c) 

Ji- 
lt was found that in some cases (Koyna, Kangra, Anjar, and Srimangal earthquakes) the 
intensity data fits a pair of straight lines instead of one, indicating thereby two rates of 
intensity decay, faster for the near distance and slower for the larger distance range. The 
reason behind this is attributed to the fact that the fault plane may be extending quite deep 
and as a consequence of which, two ends of the fault may be acting as two independent 
source of energy. There can be another possibility that the two intensity decay slopes are 
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in reality due to two independent earthquakes occurring almost simultaneously, separated 


by a few seconds, one with shallow focus and the other with a comparatively deep focus. 
2.6.6 Chandra (1979) 

Attenuation relationships were obtained using data from 31 earthquakes in and 
around the United States. As proposed by Howell and Schultz (1975), the region was 
divided into three attenuation provinces. A graphical method for the estimation of an 
initial set of epicentral intensities (/„), from the intensity-distance plot for different 

earthquake was used, thus avoiding the need to equate the maximum reported or mapped 
intensity ( ) with the epicentral distance. Using an iterative least-square fit procedure 

the attenuation relations were determined, wherein an initial approximate estimate of 
for each earthquake is successively improved. The standard equation of the type 
/ - 6 (i? + jD) -I- c log,„ {R -t- D) was used for regression. The value of D was taken 

to be equal to the average focal depth of the concerning region. The empirical relations 
between felt area (A) and epicentral intensities were also obtained. The equations derived 
were, 

San Andreas Province 

1-1^=2.065 - 0.00594 R- 2.065 log, o{R^-10) - (2.5a) 

<7=0.266 7?< 330km. 

/og,,(A) = 7.6;-h(?.72^/„ (2.5b) 

7.2 </, <10.0, cxa = 0.\9 

Cordilleran Province 

/-/, =2.819-0.00503 R-2.0J7 log, o{R + 25) (2.5c) 

<7=0.245 i?< 335km. 


log, g{A)= 1.94 + 0.448 / 


(2.5d) 
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6.8 </, <9.2, o-„ = 0.22 

Eastern Province 

/-/„ =3.374-0.00164 R-2.4141ogio{R + 25) (2.5e) 

cr= 0.363 i?< 475km. 

log,„{A) = 2.20 + 0.453I^ 

6.9 </„ <10.2, cj„ = 0.48 

Central United States 

I-l^ =3.534-0.00164 R-2.528 log, o{R + 25) (2.5g) 

C7= 0.243 R< 1600km. 

log,a{A) = 3.76 + 0.2851^ (2-5h) 

6.1</„<817, (T„ = 0.30 

2.6.7 Chandra et al. (1979) 

Data from 1 2 earthquakes from Iran was considered and equations were obtained 
for parallel and transverse to the direction of fault rupture, and for average attenuation. 
The average curve method to get the initial epicentral (/„) values, with one step iterative 

iteration method was adopted. The equations derived were, 

Along major axis 

/ - /„ = 4.8724 - 0. 00548 R-2. 708 log,o [R + 20) R<1 60 km (2.6a) 

Along minor axis 

I -1^=8.729 - 0.01158 R-6.709log,o{R + 20) i?<160kmand (2.6b) 

Average attenuation 

/-/, =6.453 -0.00121 R-4.960log,o{R + 20) i?< 120km. 


(2.6c) 
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2.6.8 Malone and Bor (1979) 

This study was aimed towards prediction of intensities at various observation 
sites, which was then compared to the actually observed intensity. The difference between 
the predicted and observed intensity was then used as site correction factors. The 
empirical relations of Evenidon (1975) were used to determine the size and depth for each 
earthquake and the local attenuation factor {k), for two physiographic parts of the United 
States. The value of k m the Puget Sound region and Northern Canada is 1.75, while 1.5 
IS more appropriate for Eastern Washington and Northern Oregon. 

2.6.9 Vanmarcke and Lai (1980) 

With data from 83 earthquakes in the Philippines, an intensity attenuation 
equation was derived in the form, 

/-/„ = 4.01 - 0.015 R - 2.10 log, o{R). (2.7) 

This compares very well with results for a tectonically similar region, the San Andreas 
region, derived by Howell and Schultz (1975). 

2.6.10 Chandra (1982a) 

Isoseismal maps for past 57 earthquakes from different parts of Greece were 
analyzed to develop relationships between intensity, magnitude and distance. The region 
under consideration was divided into six attenuation provinces, namely Hellenic, 
Peloponnese, Corinth-Saronikos, North Greece, Northwest Greece and Albania, and 
Northeast Greece province. Also, the earthquakes associated with the Hellenic province 
were grouped into two categories. Hellenic province (shallow focus), which includes 
shallow focus earthquakes showing rapid attenuation of intensity with distance, and 
Hellenic subduction zone (intermediate focus), which includes intermediate focal depth 
earthquakes, about 120 km deep, showing a rather slow attenuation of intensity with 
epicentral distance. Equation of the type / = n + i M + c /ogyg(i? + £))was 
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used. The values of the regression constants a, b, and c were found by multiple linear 


regression fit of the data. The equations derived were, 

Hellenic province (Shallow focus) 

Along Major Axis 

1 = 5.22 + 1. 154 M-0. 00394 [R + 15)- 3.3 log,, {R + 15) (2.8a) 

=0.392, (7^=0.384, cr, =0.549 for 250km. 

Along minor Axis 

1 = 5.96 + 0.890 M - 0.00565 {R + 15)- 3.3 log,, [R + 15) (2.8b) 

=0.500, (7^=0.528, (7^=0.727 for i?<165km. 

Average 

1 = 5.52 + 1.030 M - 0.00489 {R + 15)- 3.3 log„{R + 15) (2.8c) 

(7„ =0.395, =0.418, <7^=0.395, for /?<200km. 

Hillenic subduction zone (Intermediate focus) 

Parallel to the Hellenic arc 

1 = 7.11 + 1.041 M - 0.00351 R - 3.3 log „{ r) (2.8d) 

cr^, =0.697, (7;, =0.511, =0.864 for i?<460km. 

Transverse to the Hellenic arc 

1 = 6.87 + 1.213 M - 0.01222 R - 3.3 log,, (r) (2.8e) 

(7„ =0.796, fT,, =0.689, c7, =1.053 for /?< 300 km. 

Average 

1 = 6.85 + 1.137 M - 0.00713 R - 3.3 log,, {R + 15) (2.8f) 

c7„ =0.731, O'/, =0.594, (7^=0.942 for i?<375km. 


Peloponnese province 



Chapter 2 Review of Existing Attenuation Relationships 


21 


/ = 6.95 + 0.906 M - 0.00 J 369 {R + 15) - 3.3 log,, (i? + iJ) (2.8g) 

=0.226, a, =0.254, a, =0.340 for i?<80km. 

Corinth-Saronikos province 

I = 3.97 ^1374 M - 0.01157 {R + 15) -3.3 log„{R + 15) (2.8h) 

(7„ =0.465, cr, =0.615, a, =0.770 for i?<165km. 

North Greece province 

I = 5.99 + 1.014 M - 0.01037 [R + 15) - 3.3 log,, {R + 15) (2.8i) 

=0.432, a, =0.320, =0.538 for i?<160km. 

Northwest Greece and Albania province 

I = 6.29 + 0.993 M - 0.01512 {R + 15)- 3.3 log,, [R + 15) (2.8j) 

(7„ =0.362, =0.137, (7^=0.387 for i?<105kni. 

Northeast Greece province 

/ = 5.14 + 1.204 M - 0.00647 {R^- 15)- 3.3 log,, [R + 15) (2.8k) 

(7„ =0.368, a„=0.\9\, a^=0.4\5 for 7?<190km. 

This study of seismic risk in Greece, suggested that the strength of the maximum 
potential earthquakes from various seismic sources, could' be more appropriately 
expressed in terms of magnitude, rather than epicentral intensity. 

2.6.11 Chandra (1982b) 

Based on the general consideration of the regional geology and tectonics, the 
region (Indian subcontinent) under investigation was broadly divided into eight 
attenuation provinces. These are Makran, Quetta, Hindukush-Pamir, Jammu and 
Kashmir-Himachal Pradesh, Ganga basin. Peninsular province. Northeast province, and 
Bumia province. A total of 24 earthquakes were considered. The deterministic method of 
analysis was used. Equation used is of the type J-I^ = a log,,{D)+ b R-t- c log,, {R + D). 
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Where the term D is suitably assumed to be equal to the average focal depth of 
earthquakes in a particular province. The introduction of the temi D in the above equation 
will eliminate the singularity of the equation at distance i? = 0. One step regression 
analysis was performed to obtain the expression for each province. Initial set of epicentral 
intensity ( 7,^ ) values were found by the graphical method. The equations derived ai'e, 

Makran Coast 

I-I„ = 3.221 - 0.00369 R - 2.553 log,, {R + 20) (2.9a) 

cr= 0.054 i?< 600 km. 

Quetta Province 

/-/„ = 2.987 - 0.00737 R - 2.296 log„{R + 20) (2.9b) 

(7=0.429 i?<325 km. 

Hindu Kush - Pamir 

1-1,= 6.985 - 0.00590 R - 2.982 log,, [R + 220) (2.9c) 

cr= 0.301 7?<550km. 

Jammu and Kashmir - Himachal Pradesh 

1-1, = 3.975 - 0.00100 R - 3.055 log,, {R + 20) (2.9d) 

(7=0.472 7?<650km. 

Ganga Basin 

1-1, = 3.470 - 0.00210 R - 2.667 log,, {R + 20) (2.9e) 

(7-0.193 7?<1110km. 

Northeast Province 

1-1, = 2.501 - 0.00452 R - 1.922 log,, {R + 20) (2.9f) 

(7= 0.244 R< 1050 km. 


Peninsular Province 
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/-/„ = 4.987 - 0.00204 R - 3.832 log,o [R + 20) (2.9g) 

a= 0.338 i?<400km. 

Burma Province 

1-1^, =2.116- 0.001425 R - 1.626 log,,{R + 20) (2.9h) 

o-= 0.350 i?<375 km. 

It was concluded that the decay of intensity is least for the Hindukush province. This was 
attributed to the larger focal depths in this province (generally around 200-250 km). The 
peninsular and Bunna provinces show the highest attenuation rate. 

2.6.12 Rao et al. (1982) 

One-step regression analysis was performed to obtain the attenuation relations 
using data from eight Indian earthquakes. A linear relation between logarithm of intensity 
and the epicentral distance similar to that by Kaila and Sarkar (1978) was assumed. For 
each earthquake the slope of the best-fit curve was obtained. To draw the intensity 
attenuation characteristic curve for different focal depths, the relationship derived by 
Kaila and Sarkar (1978) was used. Relationship between acceleration and intensity are 
also derived. The epicentral acceleration was plotted on a semi-logarithmic scale and 
intensity on a linear scale, and the best fit straight lines obtained, were as follows, 

log (/l) = 0.233 1 - 0.05 , for firm ground, and (2. 1 Oa) 

log i„{a) = 0.288 I - 0.06 , for unconsolidated Alluvial. (2.10b) 

Where acceleration, is in the percentage of acceleration due to gravity and intensity, /, 
is on MM Scale. This study mainly resulted in the preparation of master curves involving 
the focal depth, intensity, and average peak ground acceleration of earthquake in the 
absence of strong motion data. 
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2.6.13 Benjamin (1983) 

Based on the assumption that the earthquake energy is radiated from a point 
source and subsequently transmitted through a homogeneous media, the attenuation of the 
Rossi-Forel intensity with respect to both epicentral and hypocentral distances was 
studied for the Philippine region. The data consisted of earthquakes from years 1913 to 
1980. They tried to fit many equation of the form =a + b M . The equations derived 
were, 


1 = -2.53 + 1.26 M -0.009 A, 

(2.11a) 

1 = 3.64 + 1.17 M -1.46 ln{A), 

(2.11b) 

1 = 0.91 + 1.25 M -0.86 ln{A)-0.004 A , 

(2.11c) 

I = -1.82 + 1.11 M -0.007 R, 

(2.1 Id) 

1 = 4.54 + 1.08M- 1.51 ln{R), and 

(2. lie) 

I = 1.31 + 1.13 M -0.77 ln{R)-0.004, 

(2.11f) 


The Equation 2.1 Ic is found to be the best suited for the Philippine region. 

2.6.14 Yarar era/. (1983) 

In this study fourteen earthquakes from Turkey were considered. The non- 
isotropic characteristic of the intensity distributions was accounted for, by separately, 
treating the earthquakes associated with strike slip faulting and those associated with 
Graben type. For the derivation of intensity attenuation relationships associated with 
strike-slip events, sets of intensity versus distance data points were obtained for each 
isoseismal contour, along the perpendicular and transverse axis fault rupture. One-step 
regression analysis was performed. The following results were obtained for the 
attenuation of intensity in transverse direction to the strike slip faults, 

1^-1 = 1.237 + 0.004 R + 1.21 6 ln(R) cr=0.53, and (2.12a) 


7 = 0.34 + 1.54 M - 0.001 R-1.27 ln{R) <7=0.63 


(2.12b) 
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For earthquakes associated with Graben type faulting, exhibiting more or less circular 

intensity distributions, the attenuation of intensities in radial direction from the 
macroseismic epicenter was derived as, 

I,,- 1 = 2.465 + 0.0003 R + 1.23 5 ln{R) <7=0.54, and (2.12c) 

I = 0.568 + 1.534 M -0.001 R- 1.235 ln{R) <7=064 (2.12d) 

where /„ and 1 represent, the epicentral intensity and mean intensity at a distance R in 

transverse direction to the fault respectively, and a is the standard deviation of the 
dependent variable. 

2.6.15 Burton et al. (1985) 

A set of 101 earthquakes between years 1727-1979 in Britain was considered. 
They investigated intensity attenuation with distance and focal depth, and inferred that the 

f R^ 

equation of the type /„ - 1 = a log^g — +b{R-h) fits best for the British earthquakes 

\h) 

and the equation derived was, 

( ).\ 

1,-1 = 3 log,, - +0. 0089 {R - h) . (2.13) 

\bj 

In deriving the equation one-step iterative regression analysis and average graphical 
method were adopted to get the initial 1, values. 

2.6.16 Lee and Trifunac (1985) 

Four regression equations describing the attenuation of MM intensity in tenns of 
local earthquake magnitude, M, epicentral distance, R, focal depth, H, fault size, s, and the 
local geological conditions at the recording site (-S' = 0 for alluvium, -S' = 2 for basement 
rock, and -S = 1 for intennediate sites) have been investigated. The equation used was of 
the form, 

I = 1.5 M-A-Bln{A)- — -Ds, 

^ ’ 100 


(2.14) 



Chapter 2 Review oj Existing Attenuation Relationships 


26 


where /is the MM intensity and A = {R‘ +H' + S' . For i? = 0, //=10tol5km and 
S = 0, the results of this equation become equivalent to 7 = 7.5 M -1.5 , proposed by 
Gutenberg and Richter (1956). The predicted intensity / was found to be, about one half- 
intensity level higher for alluvium sites than that of basement rock sites. 

2.6.17 Tilford and Chandra (1985) 

Data from three earthquakes of Greece was considered to obtain the intensity 
attenuation relationship. One-step iterative regression analysis was used and average 
graphical method to get the initial set of /^values was adopted. The intensity residual 

method was used to capture the effect of local geology on the observed values of 
intensity. It was observed that under similar circumstances, sites located on soil 
experiences about one intensity unit more shaking than sites located on rock. The 
equation derived were. 

Longitudinal attenuation 

/-/„ = 2.150 - 0.00750 R-2. 150 logjo{R + 10) (2.15a) 

a = 0.359 7?<300km. 

Transverse attenuation 

I-p =3.180-0.00458 R-3.180logjo{R + 10) ' (2.15b) 

a= 0.234 R<215km. 

Average attenuation 

/-/„ = 2.541 -0.00727 R- 2.541 log , o{R + lO) (2.15c) 

<7=0.233 7?<250km. 

2.6.18 Drakopoulos and Stamelou (1986) 

Intensity-distance relations along major and minor axis of proposed elliptical 

isoseismals in Western Greece considering past 41 earthquakes were developed. The 
western Greece earthquakes were divided into two attenuation provinces depending upon 
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the range of focal depths, one with focal depth in the range 0-15 km and the other in the 
range 16-41 km. The equation of the type /-/„ = a + b R + c log + D) was used. The 
values of regression constants a, b, and c were determined by least square fit between 
(/-/„) and R. Here, / is the intensity at an epicentral distance i?, D is a constant to be 
taken equal to the approximate average focal depth. The graphically estimated epicentral 
intensity set ( ) was used as first approximation. The attenuation coefficient (y) was also 

calculated using the fonnula, 

I ^ 

= ylog,a,U^ . (2.16a) 

V /r 


The intensity-distance relations under the assumption of equivalent circle and along the 
major and minor axis for earthquakes with focal depth in the range 0-15 km were as 
follows. 


/ - /„ = 3.906 - 0.004 R - 3.883 log,o {R + 10.5) 

for = 

(2.16b) 

I -I, = 1.487 - 0.013 R - 1.950 log,o {R + 5.78) 

for = 3.45 and 

(2.16c) 

1-1,^ 2.667 - 0.009 R - 3.163 log,, [R + 6.94) 

for R„„j = 3.25. 

(2.16d) 

For earthquakes with focal depths from 16-40 km: 



I - 1, = 9.289 - 0.003 R - 6.545 log,, {R + 26.3) 1 

'or /?„„„,j = 3.68. 

(2.16e) 

/ - /„ = 2.414 - 0.011 R - 1.885 log,, {R + 20.l) 

for /?„,„j = 3.96 and 

(2.16f) 

I - I, = 9.152 - 0.004 R - 6.733 log,, (R + 22.86) 

for R y = 3.53 . 

(2.16g) 


2.6.19 Gupta (1986) 

A probabilistic model, which makes use of MM intensity data, was developed for 
evaluating seismic risk for Northern India. This model predicts the probability of 
exceedance of any specified earthquake ground motion amplitude from the entire 
earthquake expected during a selected future period in the region around the site of 
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interest, and does not use more than one single design earthquake. The intensity 
attenuation model used was of the form given below, which represents the probability of 
the intensity I , being observed at an epicentral distance less than or equal to R. 


P(R) = 


loo R 


^271 a 


j e.\p , 


if pi - x^' 
a 


clx , 


(2.17a) 


where pi is the mean and cr the standard deviation of log i^R . The probability of 


observing an intensity value I j at a distance R is given by, 

P{/ = 7, } = P{/ < /, } - P{1 < [I, - 1 )] . (2. 1 7b) 

In this study 18 past Indian earthquakes were considered. The intensity attenuation 
relationship with distance were found to be as follows, 

For mean 


/ - 7, = 2.256 - 0.0099 R - 1. 798 log,o R • 


(2.17c) 


For mean plus one standard deviation 


7 - 7„ = 2.475 - 0.0048R - 2.080 log,„R . 


(2.17d) 


2.6.20 Grandori et al. (1987) 

To develop the intensity-epicentral distance relationships, data from past 19 
Italian earthquakes were considered. The epicentral distance R was considered instead of 
the hypocentral distance zl because earthquakes were of shallow focal depth range. They 
used equation of the type, 


/„ -/ = 


Iniyj) 


1 + 


1 / 7-1 






L_/ 


(2.18a) 


where 


V„ = 


Rj Rg 
Ro 


= Constant , 


(2.18b) 
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R — R 

I// = - Constant = if/ i>l. 


R-R 


(2.18c) 


/-/ 


The values found were as follows, 
ijj,, - 1.07 , y/ = 1.58 , 

where /?, is the average radius of the A isoseismal 


t^o{i„=io) — 


2.6.21 Chavez and Castro (1988) 

Two relations were proposed to predict the attenuation of Modified Mercalli 
intensity (/) with epicentral distance (D) for Mexican earthquakes. These are, 

f n \ 

ln{l) In — + 5, (D - D') + B^ ln{M A, and (2. 1 9a) 

\D' J 

f f) \ 

ln{l)=B,+B, — +BJn{D-D')+Bpn{MA, (2.19b) 

\D' J 

where M, is the surface-wave magnitude, D' is a distance related to the maximum 
intensity mapped for an earthquake. The coefficients B,, for /=0,1,2,3 were obtained by 
least-square fit, for the information contained in the intensity maps of 32 events of the 
region. The earthquakes were classified in three groups according to their epicentral 
location, focal mechanism, and focal depth. They are. 

Group 1: 18 earthqukes belonging to subduction zone, with thrust mechanisms and 
depths between 15 and 20 km. The equations derived were, 

f r)'\ 

ln{l)= 1.1090 -0.1399 In -0.0011 [D - D')+ 0.5209 ln{M A tr= Q.l\ 

\D' J 


(2.20a) 


^ D'^ 


\D'j 


ln{l)^ 1.3891 -0.0475 ~ -0.0220 ln{D- D')+ 0.3627 ln{M A cr= 0.67 


(2.20b) 
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Group (2): 10 intennediate-depth earthquakes in south-central Mexico, with normal 
faulting mechanisms and depths varying from 65 to 150 km. The equations derived were, 

f 

1.5188 -0.0627 In 0.0021 {D - D')+0.3214 ln{M^) cr=0.95 

yD'J 


(2.20c) 


f D\ 

ln{l) = 0.601 3 -0.0337 —\- 0.0224 ln{D-D') + 0.7745 ln{M^) cr=0.S0 

\D’J 


(2.20d) 


Group (3): Four earthquakes with epicenters located on the Trans-Mexican Volcanic 
belt, with reverse and nonnal mechanisms and depths shallower than about 20 km. The 
equation derived was. 


ln{l) = 3.0021 -0.3057 In ~ +0.0019 {D-D')-0.325 ln{M (7=0.83 

\D’ ) 


(2.20e) 


ln{l)= 2.0922 -0.0881 


U'J 


- 0.0233 ln{D -D')+0. 0351 ln{M ^ ) 


(7=0.79 


(2.20f) 

2,6.22 Grandori^^a/. (1988) 

It was observed that the intensity attenuation of the 
form/-/„ =a + bR + clog,„{R), does not comply very well with Italian data, as it is 

implicitly based on the assumption that energy is radiated from a point source. This 
approach is not reliable, where distance is not large compared to the source dimensions. 
Moreover, in the following formula the intensity decay (Ig-l) does not depend on the 
epicentral intensity /„, while the average trend of Italian earthquakes shows that the rate 
of attenuation is fast for small earthquakes. The drop of intensity relation used was. 
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i=I-l= 




1 + 


f-l 


Vo VK 


R 


^-1 


where 


\|/ = 


R... -R, 

V R| ~R|-i y 


R. 


Vo 


(2.21a) 


where /, is the intensity at an epicentral distance of i?, , is a constant and is 
independent of both / and epicentral intensity (I„), i?„is the minimum epicentral distance 
of the isoseismal, and R^ is the equivalent radius of the isoseismal line (the boundary 

line of the area with intensity). 

2.6.23 Gupta and Trifunac (1988) 

A probabilistic method was adopted to derive the intensity attenuation 
relationships for Indian earthquakes. The distribution of epicentral distance to a given 
isoseismal I, during an earthquake with epicentral intensity 7^, is assumed to be 
approximated by a Gaussian probability density function. The equation of the type 
/-/„ =a + b R + clog,g{R) was used. Data from 30 major earthquakes from Indian 

subcontinent were used. The equations derived for mean and mean plus one standard 
deviation for the northern India are, 

/-/„ =2.256-0.0099 R- 1.798 log, o{r), and (2.22a) 

/-/„ =2.474 - 0.0048 R- 2.080 log, o[r). (2.22b) 

2.6.24 Pietrafesa et al. (1990) 

The heterogeneous set of intensity data the equation given below is the best fit, 
while the other available equations show a good fit for Calabro-Sicillian region. 
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L. - I = 


1 




1 

+ 

1 

t 

1 

o 

1 

J] 


(2.23a) 


:th 


where and i?, represent, respectively, the equivalent radii of the first and the i 
isoseismal lines. must be evaluated with respect to epicentral intensity (/„) and the 
other terms are represented by the following ratios, 

(2.23b) 




' = Constant , and 


Rn 


Vo = 


R.., - R. 
R.-R.-i 


= Constant = \f/ 


i>I. 


(2.23c) 


2.6.25 Prasad (1991) 

It was proposed that in most of the cases, a straight line fit the intensity data points 
for Indian earthquake. This study is similar to Kaila and Sarkar (1978). The slope of these 
lines indicates decay rate of intensity with epicentral distance. The intensity-distance 
follows the relation, 

ln{l) = ln{l^)-s A, (2.24a) 


which gives ^ = exp 


-X J 


(2.24b) 


where I is the intensity at an epicentral distance of A, is the epicentral intensity. 


Intensity-distance relations were obtained for 6 earthquakes of India. These relationships 
are as follows, 

Assam earthquake (1897) 


/ 

— = exp 


-0 OOOOJ-f J 


(2.24c) 


Bihar-Nepal earthquake (1934) 


= exp 


-O.OOIOS96 


J 


(2.24d) 
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Anjar earthquake (1956) 


— -«2Xp 


-ooo.u<j:2 I 


Koyna earthquake (1967) 


= exp 


~{J OO^JcSIdi7 I 


(2.24e) 


(2.24f) 


Guwahati earthquake (1975) 





Dhramsala earthquake (1986) 


I 

— = exp 


-0 00119263 J 


(2.24g) 


(2.24h) 


It was concluded that intensity attenuation for Dharamsala earthquake is high because of 
hard subsurface structure. The intensity decay of Koyna earthquake is also high which is 
attributed to the harder nature of Deccan trap. 

2.6.26 Tao and Zheng (1991) 

The attenuation relationships used in seismic zonation of China were reviewed. 
Some recently suggested results were introduced and ongoing research using the artificial 
neural network technique to establish the relationships was reported. The seismic 
zonation in China is proposed to be guided by two principles; (1) at certain locality, the 
highest earthquake intensity experienced in the past would be also expected in the future; 
(2) regions of similar tectonic structures would be subjected to similar intensity. A total of 
566 isoseismals from 199 earthquakes were chosen as database for regression analysis 
with earthquakes having magnitude not less than 5. The attenuation relationships fitted 
for Eastern and Western parts of China were as follows. 

For eastern part of China 


T.nnp-itiidinal direction 
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I -1^=6.046 + 1.480 M-2.0811n{R + 25) a = 0.49, (2.25b) 

Transverse direction 

I -1^=2.617 + 1.425 M -1.441 ln{R + 7) a = 0.56, (2.25c) 

For western part of China 

Longitudinal direction 

1 = 5.642 + 1.538 M -2.1091 ln{R + 25) s = 0.64, (2.25d) 

Transverse direction 

I-!,, =2.941 + 1.263 M -1.494 ln{R + 7) a = 0.61 (2.25e) 

Some recently suggested results 

1. Probability models of intensity attenuation: Instead of commonly used empirical 
fonrmlae, conditional probability T*(/ = /| /^.r^jcan be adopted to describe the 

attenuation for given i,, and , where is epicentral intensity and is epicentral 
distance and written in the form of, 

P[l = i I , r")= F, (/- " I i, ,i + l)-F , (/ I i, , i), (2.27a) 

where F, is the conditional probability distribution function of for given ig and /, 
which are chosen as lognormal. 

2. A united regression analysis method: The form of the regression equation was taken 
as, 

i = A + BM -Cj lu{r„ +r^„)- C, InC + /-„* ) + e , (2.27b) 

where and are the distances for the semi-major axes and semi-minor, i.e. a and b 
mentioned above respectively; and are the predetermined coefficients for the 
two axes, and two sets of data can be obtained by the limitation of = 0 if r^^O 
and r„ = 0 if Vij ^ 0 from an isoseismal contour. 
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3. Aimed at the criticism on geological defect in the uncertainty correction procedure 


currently used in the seismic risk analysis for macroseismic intensity. The following 
form of regression equation was suggested 

y = ln{r + r,)= A + B M - C,. (2.27c) 

The intensity temi was moved to the right hand side of the equation and the tenn 
associated with distance was moved to the left-hand side as a dependent variable. It 
was done as the intensity with a decimal may not be widely accepted and VII-VI ^ 
VIII-VII, so that the sum of squares of remaining errors would be thought 
meaningless. In regression analysis the dependent variable must be continuous and it 
should be taken as part of independent variable. Thus it is more appropriate to chose 
distance as the dependent variable instead of intensity. Results showed that by this 
concept the attenuation at greater distance is faster than the conventional equations. 

4. Using the intensity-distance data, one can develop attenuation relationship by using 
the intensity-distance data; by neural network using the error back propagation 
approach, which can eliminate the error caused by smoothening of the isoseismal. 
2.6.27 Timiovska (1992) 

The data from past Macedonian earthquakes, during the years 1905 to 1985, was 
used to obtain the attenuation relation. An iterative regression analysis was used to obtain 
relation between epicentral intensity,/^, epicentral distance, R, and local geological 

conditions at the recording site, S, (<S = 0, 1 or 2 for soft, medium and hard geological site 
condition). The equation derived was of the form, 

/-/, = 0.162-0.023 R-0.577 log, o{r)-0.52S . (2.28a) 

The quality of fit was tested and the correlation coefficient computed from 


(2.28b) 
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where 1 oh is the mean of the observed intensity , and is the mean of the 

estimated intensity level 1^.^, . This study shows that the intensities recorded on sediments 

{S = 0) tend to be about one intensity level higher than the intensities recorded on the 
basement rock site {S = 2). 

2.6.28 Maugeri and Motta (1993) 

In this study 22 earthquakes from Sicily and Calabria were considered and divided 
into two groups. First group contains earthquakes with focal depth greater than or equal to 
15 km, while earthquakes with focal depths less than 15 km were kept in the second 
group. The following three type of equations were used to describe the attenuation law of 
seismic intensity. 


I = 0 h- b M + c ln(^R + Z)) 

(2.29a) 

I = a + b 1^ +cln{R + D) and 

(2.29b) 

. , r f 0,r, R 

(2.29c) 


where 1^ is the intensity of the epicentral area, R, is the equivalent radius of the highest 
mapped isoseismal line, R is the equivalent radius corresponding to intensity /, and ij / , 
are two constants obtained as follows 


Rm -Rt 

wy C 

Rt-Rt., 

for 

k>2. 

(2.29d) 

R. -R, 

'ii, 



(2.29e) 


where is the equivalent radius of the ]<!’' isoseismal line (the boundary line of the area 


with intensity). 
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2.6.29 Pantea (1994) 

Macroseismic intensity attenuation relationships for Romanian territory and 
adjacent areas were developed considering 18 earthquakes with epicentral intensities in 
the range of V to X on MSK intensity scale. The attenuation was analyzed as a function 
of distance and azimuth from the three main attenuation formulae: logarithmic, 
exponential, and power-law of which the last one was preferred, as it best fits the 
observed data. To account for the asymmetry of the isoseismal lines, the attenuation 
relationships along the following azimuths: N, NE, E, SE, SW, W and NW were sampled. 
The hypocentral distance was used where the depth of the focus was available and the 
epicentral distance in other cases. 

It was concluded that even for the zones that are adjacent to each other, there is 
noticeable difference among their rates of attenuation. Thus, future work is needed to 
verify the results obtained since it was the first approach of the macroseismic intensity 
attenuation for normal depth source on the Romanian territory. 

2.6.30 Teramoe/ a/. (1995) 

In this study an anisotropic attenuation relationship of the macroseismic intensity 
was deduced. Results obtained after analyzing a set of earthquakes in Eastern Sicily and 
Southern Calabria, showed a greater adaptability to the observed data as compared with 
those deduced using isotropic attenuation laws modified to take anisotropy into account. 

The proposed approach is based on a suitable mechanical modeling of the 
macroseismic field, which allows a quantitative characterization of its anisotropy. Taking 
into account each irregular isoseismal area and the resistance offered in the direction of 
energy propagation, it is possible to describe this resistance with the moment of inertial 
function of the generic isoseismal area , delimited by the two successive isoseismals, (i- 

if' and /'*, with reference to a coordinate system (0^), with the origin in its baricentre 
(epicenter) and with respect to the respect to the x-axis, is equal to the. 
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dA, y- , 


(2.30a) 


where A, is moment of inertia of the area , A, is the isoseismal area delimited by the 
two successive (i-1)''' and i'^‘ isoseismal, dA, is the infinitesimal element of the isoseismal 
area , y is the distance of dA^ from x-axis and ju is constant. 


The equation deduced was 


I, =I^-alog\ 




(2.30b) 


where a is the attenuation coefficient. 

2.6.31 Wang era/. (1995) 

Using the 330 intensity points from 19 historical events (Ms=3.0 to 7.8) in and 
near Jammica, a least-square regression analysis was used to develop intensity attenuation 
relationship. 

I = 2.043 + 0.99 - 1.22 log M-O.Ol R , (2.31) 


where / represents MM intensity, M ^ is the earthquake surface-wave magnitude, and R 

the epicentral distance in km. The focal depth term was not considered as the island was 
in a strike-slip plate boundary zone, so the earthquakes are almost exclusively shallow 
events (0 to 50 km). 

2.6.32 Gupta a/. (1999) 

The MM intensity data was used to study the regional attenuation characteristics 
of strong ground motion in the peninsular India, using a lognormal probability 
distribution. The best possible values of the parameters of this distribution were obtained, 
by using the isoseismal map of fifteen important earthquakes in the Peninsular India. 
One-step iterative regression analysis was performed to evaluate the attenuation 
relationship. The equation derived was. 
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/-/„ =1.8704-2.5423 R-0.001956 log, o{r), (2.32) 

where intensity (/) was measured on MM scale. It was observed that the trend of 
attenuation for the peninsular India is almost parallel to that for the Eastern United States, 
with some flat intensity-distance slope for smaller distance (horizontal translation). The 
horizontal translation of the curve for peninsular India towards smaller distance was 
attributed towards the shallower average focal depth. 

2.7 Need for the Present Study 

In India where large numbers of strong motion records are not available one has to 
rely on the isoseismals of the past earthquakes to predict the estimated ground motion at a 
particular site. The attenuation relation can be used to assess the attenuation of damage 
potential of earthquakes from epicenter. Till date very little work has been done in India. 
Only few researchers like Kaila and Sarkar (1978), Chandra (1982), Gupta (1986) and 
Gupta et al. (1999) have done some studies. These are the points, which evoke the need 
for the present study 

1 The equations fitted in the past are mostly based on the epicentral intensity (Table 
2.1). Sometimes equations involving magnitudes fit better than the epicentral intensity 
(Tao and Zheng, 1991). No studies have been made to ascertain the various equations 
based on different parameters such as magnitude and maximum intensity or both, and 
to find which is the most suitable for the Indian subcontinent. In the present study, 
equations involving magnitude, maximum intensity, epicentral intensity or their 
combination have been derived and their adequacy for the Indian subcontinent has 
been checked. 

2. In the past, equations were developed using either magnitude or epicentral intensity. 
But epicentral intensity is an imaginary intensity and magnitudes of past earthquakes 
are not much reliable as they were not measured directly from instrument. In the 
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present study, new attenuation relations have been developed considering the 
maximum intensity (intensity of innermost isoseist) and are compared with the 
conventional type of equations. 

3 Earlier equations were derived only for the average attenuation, but there exists a 
large difference in the attenuation along the longitudinal and transverse direction of 
the fault (as large as one intensity unit). Considering this point, equations have been 
developed for longitudinal as well as for transverse direction along with the average 
attenuation. 

4. The past relations describe only the mean attenuation and the deviations about the 
mean have not been taken into account. In the present study, equations have been 
developed for mean as well as mean plus one standard deviation. 

5. The past relations were developed using data points from few earthquakes only, 
which may not be reliable in all cases {e.g., Chandra, 1982(a)). In the present study 
data points from past 59 earthquakes have been considered which is expected to give 
results with more confidence. 

6. In India, equations developed in United States are still used to evaluate the seismic 
risk at a particular site. Hence there is a need to compare the attenuation of both the 
regions (Table 2.4). The past equations show that attenuation in Indian subcontinent is 
quite different from those of United States (Figure 2.2). 

7. The past equations do not directly show how damage area changes with drop of 
intensity and how the intensity drops differs from one province to other (Table 2.2). 
These things are incorporated in the present study. 


• •• 
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Table 2.1 Available equations for Indian provinces 


Province 

Researcher 

Type of equation used 

CT 

No of 
EQ. 

North 

Gupta 1986 

M, = 2 256 - 0.00990 R - 1.798 log R 

- 

18 

Chandra 1982 

I - 1, = 3 470 - 0.0021 0 R - 2.667 log R 

0.193 

one 

Northeast 

Chandra 1982 

1 - [„ = 2 501 - 0.00452 R - 1 922 log R 

0.244 

4 

Northwest 

Chandra 1982 

M, = 3 975-0.00100 R- 3.055 log R 

0.244 

3 

Peninsular 

Chandra 1982 

I - 1, = 4 987 - 0 00204 R - 3.833 log R 

0.193 

5 

Gupta 1986 

I - 1, = 1 .870 - 0.001 956 R - 2 5432 log R 


15 


Table 2.2 Available equations for United States 


Province 

Type of equation used 

a 

No of earthquake 
considered 

San Andreas 

I - 1„ = 2 065 - 0.00594 R - 2.065 log (R + lO) 

0.266 

10 

Cordilleran 

1 - !„ = 2 8 1 9 - 0 00503 R - 2.01 7 log (R + 25) 

0.245 

13 

Eastern 

I - 1„ = 3 374 - 0.00164 R - 2.414 log (R + 25) 

0.363 

8 

Central 

I - 1„ = 3 534 - 0 00 164 R - 2.528 log (R + 25) 

0.243 

10 


Table 2.3 Relative difference in expected intensity levels for geological maps units as 
shown in the geological map of California [Evernden and Thomson, 1985] 
(taken from Reiter 1991) 


Geological Map Unit 

Relative Intensity 
Compared to Granite 
Rock 

Granitic and Metamorphic Rocks 

0.0 

Paleozoic Sedimentary Rocks 

0.4 

Early Mesozoic Sedimentary Rocks 

0.8 

Cretaceous through Eocene Sedimentary 
Rocks 

1.2 

Undivided Tertiary Sedimentary Rocks 

1.3 

Oligocece through Middle Pliocene 
Sedimentary Rocks 

1.5 

Pliocene-Pleistocene Sedimentary Rocks 

2.0 

Tertiary Volcanic Rocks 

0.3 
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Table 2.4 Conversion table for intensity prepared by Kaila and Sarkar (1978) 


Oldham scale 

RF Scale 

1 

X 

2 

rx 

3 

VIII 

4 

VI-VIl 

5 

IV-V 

6 

ii-iii 

7 

I-I! 


MM Scale 

RF Scale 

I 

I 

II 

I-II 

III 

III 

IV 

IV-V 

V 

V-VI 

VI 

VI-VII 

VII 

VIII 

VIII 

VII-IX 

IX 

IX 

X 

X-XII 





Intensity' drop(M!Vl) 
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Figure 1.3 A comparison of seismic intensity scales (after Murphu and O’Brien, 1977; and 
Richter, 1958) (taken from Reiter, 1991). 
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Figure 2.4 Different definitions of distances used in the attenuation relationships [Campbell, 
1985] 



Chapter 3 

Review of Earthquake Intensity Database 


3.1 Introduction 

Isoseismal maps provide us valuable information about macro-seismic effects of 
large earthquakes. These are constrained by the field observation of the damage pattern 
and teirain changes soon after the earthquake. To develop the intensity attenuation 
relationship for Indian earthquakes, isoseismal maps from past earthquakes are collected 
from various souixes. Most of the isoseismal maps are based on the field investigations 
perfonned by the Geological Survey of India (GSI). These have been published in 
Memoirs, Records and other publications of GSI. The publications of India Meteorology 
Department such as Indian Journal of Meteorology and Geophysics, also contain 
isoseismal maps for a number of earthquakes. A few maps were also obtained from other 
sources like Bulletin of the Indian Society of Earthquake Technology and Bulletin of 
Seismological Society of America. Various workers have drawn the isoseimals using 
different intensity scales. For the purpose developing intensity attenuation relationships, 
all the maps have been converted to MM scale. 

In the present study the Indian sub-continent is divided into six provinces keeping 
in mind the general geology, tectonic setup, and the source mechanism of earthquakes. 
The six provinces are the Northwest province. Northern province, Northeastern province. 
Peninsular province, Foredeep province and Subduction province. The intensity- 
epicentral distance data for each Indian earthquake is presented in Tables 3.1 to 3.59 with 
converted intensity scale of MM. The data points considered for the present study have 
been plotted with respect to magnitude of the event in Figures 3.75 to 3.80. 
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3.2 Earthquakes Intensity-distance database 
3.2.1 Northwest Province 

This province consists of whole of northwest India, northeast of Pakistan, and the 
Kachhi plain in Pakistan. The earthquakes of this region show a mixed mechanism. The 
isoseismals of events near the Himalayas are elongated along the Main Central Thrust 
(MCT) and the Main Boundary Thrust, while those of Kachhi plains are elongated 
towards the south-southwest directions. Most of the isoseists are incomplete due to the 
veiy limited population in this region and the past isoseismals are drawn taking very few 
observation points. Some of the isoseists got accentuation due to the presence of thick 
alluvial in the Kachhi plain. The shape of isoseist are extrapolated and they are extended 
to form circles or ellipses and the corresponding distances are taken in the present study. 
The Kangra earthquake of 1905 is taken into this province because it does not fit well 
with that of north province earthquakes giving high errors. 

B. 2.1.1 Kashmir Earthquake (1885) 

This earthquake struck the Kashmir valley on May 30, 1 885 causing loss of over 
3000 lives. Jones (1985) studied the earthquake and located the epicenter near Jampur 
(34.12°N: 74.61°W), situated about 19.5 km west of Srinagar with focal depth of 12 km. 
The general trend of the isoseismals is aligned in the NW-SE direction. The epicentral 
tract covers an area of about 2641 sq. km where the maximum damage occurred (Figure 
3.1). The shock was accompanied by the sound of “a hundred canons going off at once”. 
The ai-ea over which the shock was sufficiently large to cause severe damage to buildings 
extended from the neighborhood of Srinagar on the south-east (Figure 3.2), around a little 
north of Sopur and Bramula down the Jhelum valley as far as the fort of Chikae near 
Garhi [Gosain and Arya, 1967]. 

Gosain and Arya (1967) reevaluated the damage pattern and assigned an intensity 
of IX in an area of 1902 sq. km (Figure 3.2). It appears that the maximum damage area 
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identified by Jones (Figure 3.1) encompasses the intensity IX and a part of intensity VIII. 
Jones (1885) had used the Maline scale to delineate the isoseismals, which is very rough 
and hence, this map is not used here and the isoseismals given by Gosain and Arya (1967) 
are used in the present study. The isoseismals given by Gosain and Arya (1967) are on 
MM scale but they have not explained how they have reevaluated the isoseismal map of 
.Tones (1885). 

3.2. 1.2 Kangra Earthquake (1905) 

This is one of the four great earthquakes in the Himalayan region, it took a toll of 
20,000 human lives and caused colossal loss in the form of complete damage to buildings 
and caused numerous landslides and earth fissures. The earthquake occurred on April 4, 
1905 and the epicenter was located in Himachal Pradesh (33.63°N: 76.2°E). The 
magnitude of the event was 8.6. Middlemiss (1905) studied the damage pattern on the 
basis of 1 0 point RF scale. The meizoseismal intensity, according to his interpretations, 
reached the maximum of X on RF scale (Figure 3.19) and focal depth was estimated to be 
in between 33 to 65 km. 

The elongated epicentral tracts enclosed within the three higher isoseismals were 
closer towards WNW direction indicating a rapid decrease of intensity and were widely 
separated towards ESE direction indicating a slow rate of decrease. A small isolated 
ellipse of MM intensity VII- VIII was found to exist around Mussoorie .area, which was 
attributed by Middlemiss to a separate focus. The isoseismals of the event can be obtained 
from Middlemiss (1910), Kaila and Sarkar (1978) or Dasgupta et al. (2000). Middlemiss 
prepared the basic map and Kaila and Sarkar (1978) have reproduced the same 
isoseismals intensity converted from RF scale to MM intensity scale. Dasgupta et al. 
(2000) had reproduced the isoseisnJals on the geological map of India with all geological 
and geotectonic information. 
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For the present study the isoseismals (Figure 3.19) prepared by Middlemiss (1910) 
are used because all other maps were derived from this only, and are converted to MM 
scale. 

3.2. 1.3 Baluchistan Earthquake (1909) 

This earthquake struck the Kachhi plain, alluvium basin, on October 21, 1909 
with its epicenter at 30°N: 68°E. The isoseismals are elliptical in shape with elongation 
towards southeast direction (Figure 3.3). The maximum intensity assigned was IX on RF 
scale. Elliptical shape of the isoseismal indicates that the strata beneath are uniform and 
homogeneous. The isoseismal map given by Haron (1912) has been utilized for the 
present study after converting it to MM scale. 

3. 2.1.4 North West Himalayan Earthquake (1929) 

Macroseosmic survey for this event was carried out by Coulson (1930), who 

assigned maximum intensity of VIII on RF scale, magnitude of 7.1 and computed a focal 
depth of around 200 to 250 km with its epicenter at 36.5°N: 70.5°E. The isoseismals are 
elliptical with WNW-ESE elongation (Figure 3.4). Mukherjee (1950) reevaluated the data 
of this event and constrained the isoseismals with NW-SE trend and almost unifonn 
attenuation pattern in the NE as well as SW quadrants. The isoseismals given by 
Mukherjee (1950) (Figure 3.5) have been utilized in the present study since it is directly 
on MM scale and the elongation pattern of the isoseismals are similar to the that of most 
of the isoseismals of this region. 

3. 2. 1.5 Mach Earthquake (1931) 

The earthquake occurred on August 27, 1931, felt over a vast area of 9,58,300 sq. 
km and had epicenter location at 30.2°N: 61. TE. West (1934) studied this earthquake and 
constrained the isoseismal on RP scale (Figure 3.6) and assigned a maximum intensity 
VIII on RF scale. The movement of the cmst parallel to the isoseists might have caused 
the earthquake [West, 1934]. The isoseismal map shows how the isoseismal lines extend 
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northeast in the Indus valley. This is in accordance with the general observation that the 
alluvial soil experiences larger shaking than the hard rock. The irregular shape of the 
isoseismal IV, V and VI was attributed to severe manifestation of earthquake effects in 
the alluvial tracts [West, 1934]. The isoseismal map of the event prepared by West (1934) 
has been utilized in the present study after converting it to MM scale. The maximum 
intensity assigned to it on MM scale is VIII since there are evidence of exceeding 
intensity VIII on RF scale at some places. 

3. 2.1.6 Sharigh Earthquake (1931) 

Sharigh earthquake of August 25, 1931 with epicentral location as 29.8°N: 67.7°E 
felt over a relatively small area of 960,000 sq. km and had the maximum mapped 
intensity of VIII on RF scale (Figure 3.7). West (1934) delineated the isoseismals on RF 
scale. The isoseismals show N-W elongation, which does not conform to the general 
trend in this region. The same isoseismals has been utilized for the present study with 
converted intensities on MM scale. 

3.2. 1.7 Quetta Earthquake (1935) 

Macroseismic study for this earthquake of magnitude 7.5, maximum intensity of 
IX on RF scale, with epicentral location of 29.5°N: 66.8°E were conducted by West 
(1936). The isoseismals of this event represent the similar pattern as that of Mach 
earthquake 1931, with smaller felt area of about 259,000 sq. km (Figure 3.8). Both the 
earthquakes show epicentral elongation in the NE-SW directions because of the presence 
of alluvium. For the present study isoseismals given by West (1936) are used after 
converting it on MM scale. The maximum intensity assigned to this event is X on MM 
scale because there are evidences of reaching intensity X on RF scale at some places 
[West, 1936]. An isoseist of intensity X on MM scale is drawn inside the innermost 
isoseist by approximation. The outermost isoseist is neglected since it is very irregular in 
shape and does not follow the trend of higher isoseists. 
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3.2. 1.8 Hindukush Earthquake (1937) 

This earthquake occurred in the Hindukush region on November 14, 1937 and had 
magnitude of 7.2, maximum intensity of VIII on RF scale (Figure 3.9), epicentral location 
of 37.5°N: 71.5°E, and focal depth around 200 to 250 km. Coulson (1938) studied this 
earthquake and delineated the isoseismals on RF scale. Kaila and Sarkar (1978) 
reproduced the isoseismals on MM scale. The geological map for the earthquake along 
with the isoseismal is also available [Dasgupta, 2000]. For the present study the map 
given by Coulson (1938) has been utilized after converting it on MM scale. The cause of 
earthquake in this region is reported to be due to the great dislocation crossed by 
numerous secondary faults, which separate the tertiary formations from the ancient rocks. 
Forces responsible for these faults seemed to be responsible for the earthquake also [de 
Ballore, 1940]. The isoseismals are modeled as if they represent a family of half ellipses. 
The maximum assigned intensity is VIII after converting it to MM scale, since there are 
evidences of intensity reaching higher than VIII on RF scale. The isoseist VIII is drawn 
approximately inside the isoseist VII with similar shape. 


3. 2.1. 9 Pamir Earthquake (1939) 

The earthquake occurred on November 21, 1939. Its magnitude was 6.9 and 
induced similar intensities as the two earlier earthquakes in thi's region [Mach, 1931 and 
Quetta, 1935]. Coulson (1940) who studied this event estimated the focal depth to be 200 
km and epicentral location as 36.5°N: 74°E. The maximum assigned intensity was VIII 
on RF scale (Figure 3.10). Immense dislocations at the sharp change in direction of the 
Himalayan chain, must be contributing to the instability of the region. Kaila and Sarkar 
(1978) had produced the same isoseismal on MM scale. For the present study the 
isoseismals given by West (1940) has been utilized after converting it to MM scale. 
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3.2.1.10 Cliamba Earthquake (1945) 

Kapila (1959) studied the effects of this earthquake, delimited isoseismals of 
intensities IX, VIII, and VII on MM scale and located the epicenter at 32.75°N: 76.5°E. 
The earthquake had a magnitude of 6.5 and the isoseismals roughly have a north-south 
orientation, this does not conform to the general disposition of the structural elements of 
the region (Figure 3.11). Krishnaswamy (1962) explained this occurrence with the 
presence of nearly north-south trend of possible tear faults in the area. He further 
correlated the event with the origin of the Kangra orocline. Dasgupta et al. (2000) had 
reproduced the same isoseismals on the geological map of the region. The same 
isoseismals are considered in the present study. 

3.2.1.11 Badgom Earthquake (1963) 

The Badgom earthquake of magnitude 5.5 struck the Srinagar area on September 
2, 1963. Wakhaloo (1977) studied the earthquake on the basis of damage to the structures 
and other effects like development of earth fissures, and assigned a maximum intensity of 
VIII on MM scale (Figure 3.12). The epicenter of the event was located near the village 
Kolabug (33.98°N: 74.73°E). Wakhaloo (1977) concluded that the dislocation along a 
fault with probable NNW-SSE trend, existing below Karewa sediments was the cause of 
the earthquake. The earthquake was felt over an area of approocimately 100,000 sq. km. 
On account of the shallow focus (44 km), the shaking was less intense but it was strong 
enough to result in the collapse 2000 village houses mostly of mud construction [Gosain 
and Arya, 1967]. The isoseismals of the event are also available with Gosain and Arya 
(1967). For the present study the isoseismal given by Wakhaloo (1977) has been utilized 
since it is based on the extensive field survey. 

3.2.1.12 Anantnag Earthquake (1967) 

This earthquake struck the entire district of Anantnag on February 20, 1967 and 
had a magnitude of 5.7, focal depth of around 24 km, and the epicenter at 33.60°N: 
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75.28°E. The shock was also felt in the south-eastern section of Srinagar district. The 
quake was accompanied by a deep rumbling sound (lasted for 20 seconds), which was 
described as the pealing of thunder by the local populace [Gosain and Arya, 1967]. The 
damage pattern was confined to VII isoseismal and is around the village Sop, where about 
10% old and poorly constructed houses were wholly or partially wrecked. About 458 
villages were affected and 786 houses were totally damaged in Anantnag district. 
Dasgupta et al. (2000) had produced the isoseismals given by Wakhaloo (1977) on the 
geological map of the region. The isoseismals given by Gosain and Arya (1967) have 
been utilized for the present study since they constrain more isoseismals while the map by 
Wakhaloo gives only the epicentral area and the next intensity (VII). The looping part of 
the isoseismals (Figure 3.14) on the western side has been neglected and the isoseismals 
are modeled as concentric ellipses. 

3.2.1.13 Kinnaiir earthquake (1975) 

On January 19, 1 975, the entire tribal belt of Kinnaur and Lehul Spiti was shocked 
by a strong event of magnitude around 7. It was felt over an area of about 0.25 million sq. 
km and had epicenter located at 32.45°N: 78.43°E. The isoseismals are elliptical in shape 
and oriented in the north-south direction, parallel to the Kaurik-Chango fault [Hukku et 
al, 1975]. The earthquake caused considerable loss of life and different degrees of 
damage to various constructions in the area. The most common construction in stone 
masonry suffered extensive damage. Corrugated iron sheets nailed on timber arches and 
frames forming cylindrical shell and barrack type structures suffered little damage. Large 
size boulders rolling down the hill slopes resulted in severe damage to some of the 
buildings and highways [Singh et al, 1977]. A ground fissure of about 18 km length was 
formed which appeared to be a surface expression of the fault causing the earthquake 
[Gosavi et al, 1977]. 



Chapter 3 . Review of Earthquake Intensity Database 


54 


Hukku et al. (1975) assigned a maximum intensity of X on MM scale (Figure 
3.23) on the basis of the wide and extensive ground fissures of widths up to Im; Singh et 
al. (1977) assigned a maximum intensity of IX on MM scale (Figure 3.25), as a few C- 
Type structures (RCC structures) existing in the area had not suffered damage suitable to 
assign intensity IX. Singh et al. (1977) calculated the magnitude as 6.7 and depth of focus 
as 25 km. Gosavi et al. (1977) also estimated maximum intensity IX on MM scale but 
their map contained a small isoseismal of X (Figure 3.25), which they attributed to local 
accentuation caused by loose foundation conditions. They estimated the magnitude of the 
event to be 6.2 and focal depth around 30 km. The possible sources of isoseismals are 
Hukku et al. (1975), Singh et al. (1977), Gosavi et al. (1977), Kaila and Sarkar (1978), 
Dasgupta (2000). Kaila and Sarkar (1978), and Dasgupta (2000) have reproduced the 
same isoseismal map. Out of all sources, the isoseismals of Hukku et al. (1975) of 
Geological Survey of India has been utilized for the present study. 

3.2.1.14 Dharamsala Earthquake (1978) 

This earthquake of magnitude 5 struck Himachal Pradesh (32.23°N: 76.61 °E) on 

June 14, 1978. Kumar et al. (1981) studied the earthquake and assigned a maximum 
intensity of VI on MM scale. They had interpreted that the strike slip movement along an 
N-S tear, exposed west of Dharamsala, which have been displaced could have been 
responsible for the earthquake. However Dasgupta et al. (1982), after studying the focal 
mechanism indicated a thrust type of mechanism along more or less N-S trending nodal 
planes. The isoseismals are elongated towards N-S direction and are elliptical in shape. 
The isoseismal map (Figure 3.26) given by Kumar et al. (1981) is considered in the 
present study. 

3.2.1.15 Kathua Earthquake-l(Jammu) (1980) 

Two earthquake of magnsitude 5.5 and 5.4, triggering at an interval of 12 minutes, 

shook Kathua on August 24, 1980. These were named as Bhaddu-Dudwara and Lohai- 
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Malar shocks. Krishnamurthy et al. (1980), on the basis of macroseismic survey 
estimated the maximum intensity of the Bhaddu-Dudwara event to be VIII and that of 
Lohai-Malar shock to be VII on MM scale (Figure 3.15). Isoseismals delineated thus has 
a general NW-SE elongation, which more or less conform to the regional structural 
trends. For the present analysis both the events are considered separately because both 
sets of isoseismals can be separated. 

3.2.1.16 Kathua Earthquake-2 (Jammu) (1980) 

The details of this earthquake are described in the Section 3.2.15 since both the 

events took place simultaneously. 

3.2.1.17 Dharamsala Earthquake (1986) 

This earthquake had magnitude of 5.7 and epicenter location as 32.25°N; 76.3°E. 

It struck the Kangra region on April 26, 1986. Gupta et al. (1986) studied the earthquake 
and assigned maximum intensity of VIII on MM scale (Figure 3.30). The attenuation of 
intensity VII is quicker on the southeastern side than on the northwestern side while the 
reverse is true for intensities VI and V and once again the attenuation of intensity IV was 
rapid in the southeastern direction. From these intensity decay patterns Narula and Shome 
(1992) inferred that the source fault is aligned in the N55°W: S55°E direction along the 
general trend of the isoseismals but adjustments might have taken place along the north- 
south tear located west of Dharamsala township. 

In this study the isoseismal map prepared by Narula and Shome (1992) is utilized. 
The reverse attenuation of isoseismals is neglected. The isoseismals are modeled as if 
they represent a family of concentric ellipses. 

3.2.1.18 Srinagar Earthquake (1988) 

Gupta (1988) prepared the isoseismals of the event of magnitude of 3.2. The 
epicenter of the earthquake was at 34.4°N: 75.3°E as per Indian Meteorology Department. 
The trend of the isoseismals corresponds to the Himalayan trend of NW-SE and longer 
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axis extended along the course of Jhelum River (Figure 3.16). The isoseismals of the 
event is taken from Dasgupta et al. (2000). 

3.2.1.19 Chamba Earthquake (1995) 

An earthquake of magnitude 4.5 rocked the Chamba region of Himachal Pradesh 
on March 24, 1995. Pande and Sharda (1995) studied the earthquake, assigned maximum 
intensity of VII on MSK scale, and focal depth around 10 km (Figure 3.17). The same 
isoseismal map was reproduced on the geological map by Dasgupta et al. (2000). For the 
present study the isoseismals are taken from Dasgupta et al. (2000). 

3.2.2 Northern Province 

This province consists of the whole Indo-Gangetic plain. The whole area is 
covered by very thick alluvium. The earthquakes of this region are of intermediate focal 
depth (10-30 km). The uniform and homogenous strata below give rise to elliptical 
isoseismals but sometime the alluvial deposits cause accentuation in the felt intensity. 
Preliminary analysis were carried out taking the foredeep and north Indian earthquakes 
into one category, but the results showed that the two should be separated. The final 
equations after separating the two give smaller errors and hence they are considered 
separately in the present study. 

3. 2.2.1 Delhi Earthquake (1960) 

The earthquake of magnitude of 6.0 struck Delhi (28.2°N: 77.4°E) city on August 
27, 1960, causing damage and minor casualties. Mukti Nath et al. (1969) studied the 
earthquake and assigned maximum mapped intensity of VII on MM scale (Figure 3.22). 
The depth of focus was calculated to be in between 5 to 6 km. The elliptical isoseismal 
VII indicates northeast-southwest feature parallel to the junction of Alwar quartzite and 
alluvium. The probable cause of the earthquake, as given by them, was the slippage of 
rocks along the weak regions of fault planes beneath the alluvium. As the isoseismal are 
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relationship of earthquake with the fault. The isoseismals also represent some isolated 
intensities of VI. For the present studies the isoseismal map (Figure 3.29) given by 
Srivastava et al. (1980) has been utilized. The isoseismals are modeled as if they 
represent a family of half-concentric ellipses and the isolated portions are neglected. 

3.2.2. 5 Uttarkashi Earthquake (1991) 

A moderate earthquake of magnitude 6.6 (IMD) shook Garhwal Himalayan region 
on October 20, 1991. The shock was widely felt in the adjoining northern states of 
Punjab, Harayana, Himachal Pradesh, Jammu and Kashmir, and Delhi and caused 
extensive damage in the districts of Uttarkashi, Tehri, and Chamoli. The earthquake took 
a toll of 768 human lives and caused injuries to 5066 others. The strong ground motion 
lasting for over 45 seconds caused damages ranging from hairline crack to complete 
collapse of around “one hundred thousand” houses [Narula et al., 1995]. The earthquake 
was assigned maximum intensity of IX on MSK scale. The isoseismals are elongated 
along WNW-ESE directions. The isoseists are generally elliptical in shape with some 
isolated isoseists of intensity VI and V. The direction of rupture propagation is indicated 
to be towards southeast from the widening of isoseismal in that direction. 

The isoseismals of the event are available with Narula et al. (1995), Rastogi and 
Chadha (1995), GSI (1992), and Dasgupta et al. (2000). Narula et al. (1995) delineated 
the epicentral tract and isoseismals (Figure 3.33) of higher intensities on the basis of 
comprehensive damage surveys. Rastogi and Chadha (1995) have redrawn the isoseismal 
(Figure 3.34) on the basis of data given in GSI (1992). Dasgupta et al. (2000) had 
reproduced the map given by GSI on the geological map (Figure 3.35). For the present 
study the isoseismal given by GSI (1992) has been utilized. 

3.2.2.6 Delhi Earthquake (1994) 

A mild earthquake of magnitude 4 and focal depth of 4 km struck the national 
capital on July 28, 1994. Gupta and Sharada (1996) studied the earthquake and assigned a 
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maximum mapped intensity of V on MSK scale (Figure 3.36). The isoseismals are oval in 
shape and trending in NE-SW direction. The shape of the isoseismals indicates that a NE- 
SW trending fault might be responsible for the event. 

3.2.2. 7 Bhind Earthquake (1994) 

This earthquake had a magnitude of 4.8 and shook the Bundelkhand and adjoining 
parts on September 1, 1994. The earthquake was studied and investigated by Pande et cil. 
(1995). The general trend shown by the isosiesmals is elongation in ENE-WSW (Figure 
3.37), indicating that a slip along that direction might be responsible for the event. The 
same isoseismals are used which has also been reproduced in Dasgupta et al. (2000). 

3. 2.2.8 Chamoli Earthquake (1996) 

This earthquake had a magnitude of 4.5 and struck the Chamoli district on January 
23, 1996. The isoseismals given by Singh and Joshi (1996) show general elongation 
towards N25°E-S25°W direction and had maximum mapped intensity of V on MSK scale 
(Figure 3.38). The earthquake had a focal depth around 20 km and a fault trending NE- 
SW might be responsible for it. 

3.2.2. 9 Garhwal Earthquake (1996) 

A moderate earthquake of magnitude of 5.0 rocked the central part of Garhwal 
Himalayan region on March 26, 1996. Pande and Gairola (1996) conducted the field 
survey and assigned maximum intensity of VI on MSK scale (Figure 3.39). The general 
trend shown by the isosiesmals is elongation in NW-SE direction. 

3.2.2.10 Indo-Nepal Earthquake (1997) 

This earthquake had a magnitude of 5.5, focal depth of 16 km and struck the Indo- 
Nepal border area (29.8°N: 80.5°E) on January 5, 1997. The isoseismals of the event 
(Figure 3.40) are taken from Dasgupta et al. (2000). The earthquake had maximum 
mapped intensity of VI on MSK scale. The epicentral tract was elongated in the NE-SW 
direction following the general trend of earthquakes in this area. 
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3.2.2.11 Chamoli Earthquake (1999) 

A recent earthquake of magnitude 6.3 rocked Garhwal-Kumau Himalayan region 
on March 29, 1999. The United States Geological Survey (USGS) located its epicenter at 
30.5 12°N: 79.403°E near Gopeshwar, in Chamoli District of Uttar Pradesh, and depth of 
focus as 21 km. A joint study conducted by The National Society for Earthquake 
Technology (NSET), Nepal and department of earthquake engineering, Roorkee, studied 
the earthquake and drew the isoseismal using the MSK scale. The event has a maximum 
intensity of VIIE on MSK scale (Figure 3.41). The elongation of isoseismal line of 
intensity VII in WSW-ENE direction depicts that the strike of the slip may be in WSW- 
ENE direction. The elongation and anti-clockwise rotation of the isoseismal line for 
intensity VI may be due to the effect of the Alaknanda river valley in the almost NE-SW 
direction. The pattern of damage i.e., maximum damage in the strike direction, at some 
distance from the epicenter and not near the epicenter depicts that the slip has been on a 
low angle fault. The isoseismals of the event is taken from Shanker and Narula (1999). 
3.2.3 Foredeep Earthquakes 

This province consists of the Indo-Gangetic foredeep. The earthquakes in this 
region show quite different characteristic as those of Northern India province. The 
earthquakes of this province have focal depth generally larger then that of Northern India 
province. The isoseismal of this province are generally accentuated in the southern 
direction due to the presence of thick alluvial below and the northern portion are 
generally incomplete due to lack of population and difficult terrain. 

3. 2.3.1 Kathmandu Earthquake (1833) 

The epicenter of the event lies close to Kathmandu. The magnitude of the event 
was reported to be in between 7.5 and 7.9 [Dasgupta el al, 2000]. The earthquake struck 
at 1 1 A.M. by a loud subterranean noise comparable to the noise produced by a discharge 
of 100 pieces of artillery. The direction of the motion was variously stated. At Tirhoot it 
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was said to have been from east to west; at Buxer, from north to south [Oldham, 1883]; at 
Patna apparently east and west; at Calcutta, northeast to southwest; and at Kathmandu, 
apparently east to west. All shocks came from east to northeast. Each of the shock, lasted 
three or four seconds but some are stated to have lasted one minute [Oldham, 1883]. The 
event affected around a million square km region of north India, Nepal and Tibet and 
caused landslides, rockfalls and destruction of 4600 dwellings. It took around 500 lives. It 
has maximum mapped intensity of IX on MM scale (Figure. 3. 18). The isoseismals are 
elliptical in shape with the elongation towards NW-SE directions. The fall of intensity is 
faster towards the Himalayas but it was slower towards SW direction. 

3.23.2 Bihar-Nepal earthquake (1934) 

On January 15, 1934, north Bihar and Nepal were shaken by one of the severest 
earthquake in the history measuring 8.3 on the Richter scale. Its epicenter was located 
near Madhubani (26.3°N: 86.3°E), and maximum intensity observed was X on MM scale 
(Figure 3.20). The earthquake was studied by GSI (1939). The intensity was found to 
diminish rapidly with increase in the distance away from epicenter lying on the eastern 
edge of the region covered by isoseismal X. The high intensity region of the Nepal valley 
was attributed to the presence of unconsolidated alluvium. The elliptical shape of the 
isoseismals was supposed to be due to elongated nature of the fault rupture. The shock 
was felt over an area of 4,920,000 sq. km. Roy (1939) assigned depth of focus from 
seismometric study as 14.8 km while Richter (1958) kept the focus between 20-30 km 
depth. 

3.2.3. 3 Bulandshahar Earthquake (1956) 

An earthquake of magnitude 6.5 struck the Bulandshahar and Khurja towns 

(28.2°N: 77.8°E) on October 10, 1956. The isoseismal map was taken from Tandon 
(1975). The depth of focus was calculated to be 20 km and radius of perceptibility 300 
km. It was indicated that a strike slip motion along NNE-SSW was responsible for the 
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event. The isoseismals are elongated along southeast direction. The same isoseismals are 
considered in the present study. The earthquake was felt throughout northwest Uttar 
Pradesh and adjoining area and caused minor property damage at Bulandshahar and 
Khurja. It was strongly felt in Delhi where a large number of people rushed out in open 
for safety. Minor cracks also appeared in a few buildings. In Delhi the intensity was close 
to VI on MM scale (Figure 3.21). 

3. 2.3.4 Roorkee Earthquake (1975) 

A moderate earthquake of magnitude 4.7 shook the area around Roorkee on 
November 6, 1975, with the epicenter at 29.85°N: 77.88°E. The earthquake was felt at 
Delhi, Dehradun, Bijnor, and other adjoining districts in the state of Uttar Pradesh. 
Different degrees of damage were caused mostly to the older constructions around 
Roorkee. Cracks along mortar joints and fall of plaster were commonly observed in a 
number of residential buildings. The earthquake was recorded on Multiple Structural 
Response Recorder installed at the premises of the School of Research and Training in 
Earthquake Engineering, University of Roorkee. This was for the first time that an 
earthquake had been recorded on a strong motion instrument installed in Roorkee [Arya 
etal., 1977]. 

The area around Roorkee has a 300m thick cover of recent alluvium. The region 
north of the area consists of Terai and Bhabar formation, which are unconsolidated young 
deposits. Further to the north, the Indo-Gangetic plain is terminated by Quaternary 
tectonic movements. The geomorphologic evidences indicate the existence of two 
important faults near Roorkee in the unconsolidated young deposits namely, the 
Ganeshpur fault and Roorkee fault. The Roorkee fault runs just north of the terrace bluff 
of Roorkee. The apparent termination of the tilted terraces, sudden swing of regional 
drainage from N-S to W-E, apparent relative young age of the sediments in the zone north 
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of the fault and the terraces bluff are evidences of the existence of this fault which strikes 
nearly NW-SE direction [Arya et al, 1977]. 

Arya et al. (1977) studied the earthquake and assigned a maximum mapped 
intensity of VI on MM scale. The elliptical isoseismals as given by Arya et al. (1977) are 
elongated towards NW-SE (Figure 3.27), which is parallel to the Roorkee fault. Kaila and 
Sarkar (1978) had reproduced the same isoseismals. The isoseismal given by Arya et al. 
(1977) has been utilized in the present study. 

3.23.5 Biliar-Nepal Earthquake (1988) 

This earthquake struck on August 21, 1988 (epicenter at 26.72°N: 86.63°E) and 
had magnitude of 6.4 (USGS). The earthquake severely rocked the entire north Bihar and 
adjacent areas of east Nepal, Sikkim and Darjeeling in West Bengal. The shock was felt 
almost all over north and eastern Indian from Rajasthan in the west to north eastern States 
in the east and from Himalayan foothills in the north to Chhotanagpur plateau in the 
south. Darbhanga and Madhubani were the worst affected districts. The earthquake killed 
about 300 people in north Bihar, 700 people in Nepal, and caused collapse of about 
50,0000 houses in the entire area. The reason for the damage being so extensive was that 
the quake struck close to thickly populated areas. Moreover, the alluvial soil of the Indo- 
Gangetic plains makes the region particularly vulnerable during earthquakes. The quake 
lasted for 40 seconds, and its impact was a catastrophic release of energy that struck the 
surface in ever-expanding concentric circles, leaving in its wake a trail of death and 
devastation [Duggal and Sato, 1989]. 

Sinha (1993) studied the earthquake and assigned a maximum intensity of IX on 
MM scale (Figure 3.31). The isoseismals contain an isolated isoseist of intensity VIII 
towards the south of Ganga river. The mechanism solution of the event by Banghar 
(1991) indicates two nodal planes, one striking N58°W: N88°E, and the other with 
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N30°E: N60°W. The isoseismal map prepared by Thakkar et al. (1993) is symmetrical 
(Figure 3.32) and quite different from that of Sinha (1993). They assigned maximum 
intensity of VIII on MM scale. For the present study the isoseismal map given by Sinha 
(1991) has been utilized since it looks quite realistic and reliable since alluvial soil may 
cause local failure, which distort the shape of isoseismal from general elliptical shape. 
Kumar (1990) has also given a preliminary isoseismal map. The isoseismals are modeled 
as if they represent a family of half-concentric circles and the isolated isoseists are 
neglected because of local effects encountered. 

3.2.4 Northeastern Province 

The earthquakes of this province are of intermediate focal depth and the , 
isoseismals are generally elliptical in shape. The soil beneath is thick alluvial of 
Bhannputra river plain. The isoseists sometimes get accentuated due the presence of thick 
alluvium. Preliminary survey were conducted keeping the northeast and the subduction 
earthquakes together but the results showed that the two should be separated. The final 
equations after separating the two give lesser errors and hence they are considered 
separately in the present study. 

3. 2.4.1 Assam Earthquake (1897) 

The Assam earthquake of 1897 comes amongst the greatest earthquakes of the 
world. It struck Shilong (26°N; 91°W) on June 12, 1897 and had a magnitude of 8.7 and 
epicentral intensity of XII on MM scale (Figure 3.42). Oldham (1899) determined the 
focal depth in between 180 and 320 km. Oldham (1899) attributed cause of this 
earthquake to a movement along thrust plane(s) and along secondary thmst and fault 
planes. He observed dislocation up to 10.66m along Chedrang fault and only 3 cm along 
comparatively minor Samin fault. The elliptical shape of the isoseismal, given by 
Oldham (1899) were closely spaced in the east and southeast directions but widely spaced 
in the west and southwest directions. The approximately elliptical shapes suggest that as 
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if the soil underneath is homogeneous. The inner irregular isoseismals were due to the 
geological conditions, the alluvial causing a higher degree of damage. Oldham (1899) 
estimated the felt area of the event to be not less than 3,120,000 square km. 

The isoseismals of the events can be obtained from many sources like Oldham 
(1899), Kaila and Sarkar (1978), Narula and Sharda (1997), and Dasgupta et al. (2000). 
Oldham (1899) first prepared the isoseismal map of the event on the basis of field 
observations on Oldham intensity scale. Kaila and Sarkar (1978) had reproduced the map 
(Figure 3.43) with converted scale of intensity from Oldham to MM intensity scale. 
Narula and Sharda (1997) reviewed the field data available [Oldham, 1899] and 
regenerated the isoseismals on MSK scale and assigned maximum intensity of XI on 
MSK scale (Figure 3.43). Dasgupta et al. (2000) reproduced isoseismal map (Oldham, 
1899) on the geological map of the region with all geological and geotectonic 
mfomiation. For the present study the isoseismal prepared by Narula and Sarkar (1978) 
has been utilized since it is directly on MM scale. 

3. 2. 4. 2 Cal cutta Earth quake (1906) 

Middlemiss (1907) delineated the isoseismals (Figure 3.44) and located the 
epicenter at 27°N; 91°E. There is no record of estimate of magnitude of the earthquake. 
The epicentral area was elongated NS and coincided roughly with the course of Bhagirithi 
River down to its junction with Hooghly River. Kaila and Sarkar (1978) have reproduced 
the same isoseismal map on MM scale. 

3. 2. 4.3 Shrimangal Earthquake (1918) 

This earthquake had a magnitude of 7.6 and epicenter at 24.5°N; 91°E. The 
earthquake was studied in detail by Stuart and used Oldham intensity scale to draw the 
isoseismals. It was suggested that the cause of earthquake was the subsidence along the 
southern side of a fault running below the Sylhet valley alluvium, approximately under 
the major axis of the epicentral tract. The calculated focal depth was around 13 to 14.5 
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km. The isoseismals were found to be very regular with the apex towards northwest. The 
isoseismals (Figure 3.45) of the event are available with Stuart (1926), Kaila and Sarkar 
(1978), and Dasgupta (2000). For the present study the isoseismal map given by Stuart 
(1926) is used because it is derived from the actual field observations. Kaila and Sai’kar 
(1978) have reproduced the same isoseismals with converted intensity scale to MM scale, 
while Dasgupta (2000) had produced the same isoseismals on the geological map of the 
region. For the present study the isoseismals given by Staurt (1918) has been utilized after 
converting to MM scale. 

3. 2. 4.4 Dliiibri Earthquake (1930) 

This earthquake had a magnitude of 7.1 and epicentral location at 25.5°N: 90.9°E. 
It was studied in detail by Gee (1934). The maximum mapped intensity was IX on MM 
scale (Figure 3.46). The intensity of the event decreases rapidly northwards across the 
strike of the Himalayan. Gee suggested a considerable depth of focus, as there was no 
rapid decrease in intensity in the vicinity of the epicentral area. Kaila and Sarkar (1978) 
also drew the isoseismal map on MM scale and the same is utilized in the present study. 

3. 2.4.5 Cachar Earthquake (1984) 

This earthquake struck the northeast India on December 31, 1984. It had a 
magnitude of 5.5 and focal depth of 2 km. The isoseismals are taken from Dasgupta et al. 
(2000). The isoseismals are elliptical in shape. The epicentral tract is elongated in NE-SW 
direction following the general tectonic trend in the region. The epicentral tract was 
assigned a maximum mapped intensity of VII on MM scale (Figure 3.51). Mud houses 
were badly damaged and their collapse caused 13 deaths. The Assam type construction 
(made with bamboos for earthquake resistant) using brick masonry had mainly suffered 
damage to masonry portions. The light timber construction had generally tilted due to 
absence of diagonal bracing and not burying the timber posts into the ground [Agarwal, 



Chapter 3 Review of Earthquake Intensity Database 


67 


3.2.5 Subduction Earthquakes 

Subduction earthquakes have high focal depth (60 to 100 km). The high focal 
depth is mainly attributed to the subduction of Indian plate with that of Burma plate. The 
earthquakes near the Dauki and Halflong fault with high focal depth are considered in this 
province. The Assam earthquake of 1975 is also considered in this province because it 
has large focal depth and its epicenter is near the Dauki and Halflong fault and hence is 
considered in this province, but past literature does not mention weather it is subduction 
earthquake or not. 

3.2.5. 1 Assam Earthquake (1950) 

This earthquake struck on August 15, 1950 and the epicenter was near the 
northeastern border of India (28.5°N: 96.5°E). The earthquake was studied by many 
researchers. The reports are available in a Central Board of Geophysics publication 
entitled “A Compilation of Papers on Assam Earthquake”. The magnitude of the shock 
was 8.7 and focal depth of 30 km (Poddar, 1953). Ray (1953) reported that the maximum 
intensity at the epicenter might have reached the value XII (Figure 3.48). The isoseismals 
elongated towards southwest might indicate the presence of a causative fault trending NE- 
SW. The isoseismals of the event are also available with Dasgupta et. al. (2000) and the 
same are utilized in the present study since they are based on the extensive field 
investigations conducted by GSI. 

3. 2.5. 2 Assam Earthquake (1975) 

This earthquake was studied by Gosavi et al. (1977). It took place on July 8, 1975, 
had a magnitude of 6.7, epicenter at 21.5°N; 94.7°E and focal depth at about 60 km. The 
earthquake was felt over an area of 500 km radius. The maximum mapped intensity 
assigned to it was VII on MM scale. The isoseismals are elliptical in shape with epicentral 
elongation towards NE-SW direction (Figure 3.49), as shown by other earthquake in this 
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region. The isoseismal given by Gosavi et al. (1977) has been utilized for the present 
study. 

3. 2.5. 3 Nagaland Earthquake (1970) 

On July 29, 1970, north eastern part of India, near India-Burma border was struck 
by an earthquake of magnitude 6.4, with a focal depth of 59 km. the epicentre was located 
at 26°N: 95.4°E. Guha and Gosavi (1974) prepared the isoseismal map and the same is 
considered for analysis in the present study. The maximum intensity assigned was VII on 
MM scale to area around Margherita and Digboi (Figure 3.50), where some damage to 
housing property occurred. The epicentral tract runs NE-SW and indicates possible 
tectonic adjustment along Haflong-Dissang thrust fault. Distribution of earthquake 
intensities show the influence or regional geology and local foundation conditions such as 
at Imphal where local accentuation in ground vibrations could be attributed to thick 
sedimentary cover [Guha and Gosavi, 1974]. 

3.2. 5.4 Burma-India Earthquake (1988) 

This earthquake had a magnitude of 6.8 and rocked the northeastern states on 

August 6, 1988. Kumar (1992) on the basis of field investigation prepared the isoseismal 
map and assigned a maximum mapped intensity of VII on MM scale (Figure 3.52). The 
isoseismals are elliptical in shape. The elliptical shape of the isoseismal with its major 
axis having NW-SE trend do not show a close association with the existing major 
lineaments in the region. Berlongfer and Diphu villages had highest and second highest 
peak ground accelerations of 0.337g and 0.331g respectively. It was reported that the 
earthquake was also felt at distant places like parts of West Bengal, and Bihar [Kumar, 
1992]. The earthquake was felt strongly in northeast. Buildings in Guwahati, Jorhat, 
Sibsagar, Diphu, Silchar in Assam and other parts of northeastern states also suffered 
damage. The earthquake caused landslides at several places including National Highway 
37 and Lumding-Diphu section of meter gauge railroad in Assam, 
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3.2.6 Peninsular Earthquakes 

3.2.6. 1 Kutch Earthquake (1819) 

This earthquake struck the Kutch region on June 16,1819. This earthquake had 
magnitude of 8.3 [Jain, 1998]. It was one of the largest intra-plate earthquakes to have 
occurred in the world. The shock was felt over a radius of 1600 km and was perceptible 
all over the country as far as Calcutta and Madras. It provided the earliest clear and 
circumstantially occiuTence of faulting. It caused fault scrap of about 90 km long and 
about 6-9 m high, trending roughly E-W and was later called as "Allah Bund (wall of 
God)" [Jain, 1998 and Dasgupta et al, 2000], and caused a tsunami from Arabian Sea to 
surged across the Rann. The dislocation model based on deformation data suggests 
rupture along a reverse fault. The location of the epicenter of this earthquake is equivocal. 
According to many researchers the epicenter lies in the south of the Allali Bund, while a 
few consider it to be lie north of it [Dasgupta et al, 2000]. The isoseismal map is taken 
from Dasgupta et al (2000). The maximum mapped intensity was IX on MM scale 
(Figure 3.53). Isoseists corresponding to intensity IX, VIII and II only were shown in the 
map. The isoseist II is very close to other isoseismal on western side and very far towards 
southeastern side. The epicentral area was very large as compared to that of other 
earthquakes in this region. Due to these reasons this event is not considered in the present 
study. 

3. 2. 6.2 Coimbatore Earthquake (1900) 

This is one of the oldest peninsular shield earthquakes. It occurred on February 8, 
1900 and epicenter was located at 10.75°N: 76.75°E. It was studied by Basil (1964) and 
assigned a maximum intensity of VII on MM scale (Figure 3.54). The depth of the focus 
was around 70 km and the magnitude was 6. From the isoseismal pattern, it was 
concluded that a fault striking north-northwest on the western coast was responsible for 
the earthquake. The isoseismal map also contains an isolated isoseist of intensity VI. For 
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the present study the isolated isoseist is neglected assuming the existence of particularly 
localized weak strata. Kaila and Sarkar (1978) and Dasgupta (2000) have reproduced the 
same isoseismals. 

3.2.6.3 Satpura Earthquake (1938) 

This earthquake struck the Satpura region on July 14, 1938. It had a magnitude of 
around 6 and epicentral location of 21.53°N: 75.83°E. The earthquake was studied by 
Mukherjee (1942). The focal depth was determined to be 43 km and maximum intensity 
was VII on MM scale (Figure 3.55). The epicentral tract is elongated along NE-SW 
direction as if fault plane movement has occurred in that direction. For the present study 
isoseismals are taken from Dasgupta et al. (2000). 

3. 2.6.4 Anjar Earthquake (1956) 

On July 21, 1956 earthquake of magnitude 7.0 and depth of focus around 13 to 18 
km, struck near Anjar town (23.3°N; 70°E) of the Kutch region. The maximum mapped 
intensity was IX on MM scale. The earthquake was felt as far away as Bombay in the 
south, Hyderabad and Thar Parkar districts of Pakistan in the north, Lorwada railway 
station in east, and the extreme western end of Kutch in the west. The area of maximum 
damage was confined to about 2,000 sq. km in the central mainland of Kutch. The 
elliptical isoseismals given by Tandon (1959) have the major axis in the northeast- 
southwest direction (Figure 3.56), indicating that the probable fault was in this direction. 
Tandon (1959) observed that the damage to the eastern portion of the Anjar town was 
more severe than to the western part. He attributed this to the fact that the eastern part of 
the town was built on soft ground whereas the western part was on hard trap rocks. This 
earthquake caused 115 fatalities and hundreds were injured. About 1350 houses were 
destroyed in Anjar town alone and almost 200 houses were partially damaged. The quake 
also resulted in landslides, which dislocated railway services and snapped communication 
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and telegraph lines. The earthquake had radius (Figure 3.57) of perceptibility around 300 
km [Dasgupta et ciL, 2000]. 

3.2. 6.5 Koyna Earthquake (1967) 

This was one of the most destructive earthquakes in the peninsular region and 
took place on December 10, 1967. Its epicenter was located at 17.2°N: 73.72°E. Chatterji 
et cil. (1969) have given the isoseismal map (Figure 3.59), prepared from the field 
investigations conducted by GSI. Central Water and Power Research Station (CWPRS), 
India has also prepared an isoseismal map. The maximum mapped intensity was IX on 
MM scale. The pattern of the isoseismals is elliptical and aligned nearly north south 
which probably indicated the strike of the causative fault. The magnitude of the event 
ranges from 6 to 7.5 given by different agencies and focal depth around 10 km. 

3. 2.6. 6 Bhadrachlam Earthquake (1969) 

The earthquake took place on the peninsular shield area on April 13, 1969. 

Mukherjee (1971) studied the earthquake, determined the focal depth around 40 km and 
magnitude as 6.5. The epicenter of the event is located at 17.9°N: 80.6°E. The elliptical 
isoseismals were having their major axis trending almost east-west indicating that the 
causative fault may be striking in that direction. The isoseismal map is also available with 
Raju et al. (1975), Kaila and Sarkar (1978), Varma et al. (19'70) and Dasgupta (2000) 
(Figure 3.61). The isoseismal (Figure 3.60) prepared by Mukheijee (1971) has been used 
in this study. 

3. 2.6.7 Broach Earthquake (1970) 

The earthquake took place on March 23, 1970, had a magnitude of around 6 and 
maximum mapped intensity of VII on MM scale. The epicenter of the event was located 
at 21.68°N: 73.2°E. The isoseismal map (Figure 3.62) was given by Chaudhury et al. 
(1972) and reported depth of focus was in between 15 to 25 km. The same isoseismals are 
considered in the present study. The isoseismals are almost elliptical in shape with their 
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elongation in the NE-SW direction indicating as if a fault movement had taken place in 
that direction. The isoseismals of this event are also available with Chaudhury et al 
(1972) and Dasgupta et al. (2000). A ground fissure of about 15 km in length in alluvial 
tract was fomied along the southern bank of the river Narmada. Spouting of fresh and 
cool water fissured through ground, very loud sounds in the epicentral zone, disturbance 
in the Nannada river water, extreme rotatory motion experienced by residents of Broach 
were noticed in this earthquake. Thick alluvial bed of the northern bank of Narmada had 
considerably accentuated the seismic movements and contributed to these extreme local 
phenomena at the epicenter although the earthquake was of magnitude 6 only [Varma et 
al, 1970]. 

3. 2.6.8 Koyna Earthquake (1973) 

Gosavi et al. (1977) studied this earthquake at Koyna, it took place on October 17, 
1973. It had magnitude of 5.2 and maximum intensity of VI on MM scale. It was felt over 
an area of 1 50 km radius and had a focal depth of 10 km. The epicenter of the event was 
located at 17.7°N: 73.9°E. The isoseismal pattern is almost elliptical in shape and 
elongated in NE-SW direction. For the present studies the isoseismal map (Figure 3.64) 
given by Gosavi et al. (1977) is considered and the outer isoseist are extended to make 
full elliptical shape. 

3. 2.6.9 Shimoga Earthquake (1975) 

There were no past records of earthquake activity in region surrounding Shimoga 
before the Shimoga earthquake of May 12, 1975. The shock was studied by Gosavi et al. 
(1977). They assigned a magnitude of 5.0 and maximum mapped intensity of V on MM 
scale (Figure 3.65). The epicenter was located at 13.8°N: 75.3°E. The isoseismals were 
found to of crescent shape with the concave side facing north-northeast. Focal depth and 
felt radius as reported by them were 35 km and 140 km respectively. 
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3.2.6.10 Jaisalmair Earthquake (1991) 

This earthquake struck the Rajasthan state on September 8, 1991 and the epicenter 
was located at 26.4°N: 70.66°E. It had a magnitude of 5.5 and focal depth around 30 km. 
The isoseismals (Figure 3.66) used in this study are taken from Dasgupta et al. (2000). 
The isoseimals given are elliptical in shape with elongation towards NW-SE directions. 
The intensity falls rapidly from intensity level VIII to VII but slowly from VII to V. The 
maximum mapped intensity assigned was VIII on MM scale. 

t 

3.2.6.11 Saurastra Earthquake (1993) 

This earthquake struck the Saurastra coast on August 24, 1993. The epicenter was 
located at 20.7°N: 71.44°E with focal depth of 29 km. The isoseismals (Figure 3.68) used 
in this study are taken from Dasgupta et al. (2000). The epicentral tract is elongated 
towards NE-SW direction but the lesser isoseists are elongated towards E-W direction. 
The maximum mapped intensity was V on MM scale. 

3.2.6.12 Killari (Latur) Earthquake (1993) 

A devastating earthquake of magnitude 6.4 (USGS) occurred in the early morning 

of September 30, 1993 in the southern part of the Latur District of Maharastra State. It 
had epicenter at 18.07°N: 76.62°E, and focal depth around 12 km (IMD). The earthquake 
caused widespread damage in some 80 villages and took a toll of more than 1 1,000 lives. 
Historically, no earthquake is reported to have occurred in this area [Gupta, 1994]. With 
death toll exceeding 8933 human lives, the Latur earthquake is globally the deadliest 
earthquake in the stable continental crust to have occurred in the recorded history. The 
focal Mechanism solution obtained by NGRI with first motion directions from regional 
and teleseismic stations indicates reverse faulting with NW-SE orientation of the nodal 
planes [Ramakrishnan et al, 1995]. 

Ramakrishnan et al (1995) has drawn the isoseismals on the basis of extensive 
field investigation conducted by Geological Survey of India. The maximum mapped 
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intensity assigned was IX on MSK scale. The higher isoseists are closer towards NE and 
wider on SW direction. Dasgupta et al. (2000), Figure 3.68, and Narula et al. (1996), 
Figure 3.69, have also given the isoseismal map. Sinvhal et al. (1994) had given 
preliminary isoseismals of the earthquake. For the present study the isoseismals given by 
Ramaki'ishnan et al (1995) are used since it is based on the of extensive field 
investigation conducted by Geological Survey of India. 

3.2.6.13 Bonaigar Earthquake (1995) 

This earthquake had a magnitude of 4.4, focal depth around 21 km and epicenter 
at 21.66°N: 85.59°E. The isoseismals (Figure 3.70) used in this study are taken from 
Dasgupta et al. (2000). The maximum mapped intensity was V on MSK scale. The 
epicentral tract is elongated towards NE-SW directions. The intensity falls rapidly in NE 
direction but slowly in SW direction. 

3.2.6.14 Jabalpur Earthquake (1997) 

This earthquake struck the Jabalpur city on May 22, 1997. Its epicenter was 
located at 23.083°N: 80.08 1°E (USGS). The focus was at a depth of 36 km. This 
earthquake was first to hit any populous city in India. This earthquake caused 
considerable damage and loss of life. Jabalpur city and several villages located within a 
radius of about 15 km around the epicenter suffered heavy damage [Rai et al., 1997]. 
Jabalpur is situated in Narmada rift valley, which is bounded, by Narmada north and 
Narmada south faults passes through the southern part of Jabalpur city along the Narmada 
river course. The deep seated earth settlement (37 to 43 km) below the surface for crust 
part of earth along the fault zone, has been considered as the reason for Jabalpur 
earthquake. 

Rai et al. (1997) studied the earthquake and assigned maximum intensity of VIII 
on MM scale (Figure 3.71). The epicentral tract was elongated in NEN-SWS direction. 
The shape of the isoseismal depicts shallow nature of the earthquake. The change in the 
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orientation of the major axis of elliptical isoseismal lines from NEN-SWS to ENE-WSW 
with the decrease in intensity can be attributed to the effect of local geology and rupture 
direction [Rai et ah, 1997], Chaturvedi et al. (2000) obtained the isoseismal map (Figure 
3.73) based on the extensive field survey conducted by GSI. The isoseismals show NE- 
SW elongation with maximum intensity of VIII on MSK scale. For the present study the 
isoseismal map prepared by Chaturvedi et al (2000) is used since it is based on the 
extensive field survey. 

3.2.6.15 Bhuj earthquake (2001) 

A devastating earthquake of magnitude 7.7 struck Kutch region of Gujarat on 
January 26, 2001. The earthquake has taken huge toll of life and property in the 
meizoseismal area as well as in cities of Ahmedabad, Bhavnagar, and Surat located at 
large distances from the epicenter. The maximum intensity was X on MSK scale. The 
epicentral tract was elliptical with an area of about 2100 sq. km, which was characterized 
by adobe, brick and stone masonry buildings and many reinforced cement concrete 
buildings. The epicentral region had suffered damage of grade 4 to 5 (MSK damage 
grades). Wide ground fissures and collapses of low height road embankments were 
observed. The long axis of the epicentral tract was aligned in ENE-WSW to E-W 
direction. 

This earthquake was first to hit major metropolitan cities and the modem 
industrial constmction of the country. In Ahemadabad, 69 RCC buildings collapsed while 
cities of Rapar, Anjar, Bachho and Gandhidham were completely destroyed. Widespread 
liquefaction was observed in the Rann of Kutch, the Little Rann of Kutch as well as the 
coastal areas in the vicinity of Gandhidham, Kandla, and between Malya and Samakhiali. 
The liquefaction was recorded profoundly in isoseist VIII to X. The distribution deviated 
from the general E-W elongation consistent with higher isoseismal, to almost north-south 
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patter along the area occupied by deep seated Quaternary and Cenozoic cover sediments 
along the Cambay Graben. The accentuation of the motions because of these thick covers 
has resulted in modification of boundaries (Figure 3.75) of isoseist VII [Narula and 
Chaubey, 2001]. The isoseismals of the event are available with Narula and Chaubey 
(2001) and GSI (2001), (Figure 3.76). For the presents study the isoseismal given by 
Narula and Chaubey (2001) are used. There are modifications in the depositions and trend 
of VII isoseist of this earthquake because of the presence of thick alluvial in the Cambay 
Graben, which has resulted in intensity of more than VII in Ahemadabad. For working 
out the attenuation the average radius of this has been considered in which the effect of 
accentuation is not there. 


• •• 
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Earthquakes Intensity Data 

Symbols used 

In. : Intensity in MM scale 
Rad. : Equivalent radius of circle 

Long. : Distance in longitudinal direction of isoseismal from geometric center 
Tran. : Distance in transverse direction of isoseismal from geometric center 
So. : Source of the isoseismal 

North-West Earthquakes Intensity Data 


Table 3.1 Kashmir earthquake (1885) 


Kashmir- 
1885 (M7) 


So. 

Jones (1885) 

Gosain and Arya (1967) 

In. 

Area 

Rad 

Long 

Tran 

Area 

Rad 

wmtm 

Tran 

IX 

2641 

29 

45 

20 

1902 

25 

63 

39 

VIII 

13678 

66 

65 

73 

4070 

36 

92 

57 

VII 

284487 

301 

363 

235 

- 


- 

- 


Table 3.2 Kangra earthquake (1905) 


1 

So. 

Ml 

iddlemiss 1910 

In. 

Area 

Rad 

Long 

Tran 

XI 

520 

13 

10 

10 

VIII.V 

4659 

40 

46 

25 

VII 

12868 

64 

74 

40 

v.v 

548912 

418 

480 

365 

IV 

1368478 

660 

725 

600 

II.V 

4154756 

1150 

1250 

950 


Table 3.3 Baluchistan earthquake (1909) 



So. 

Haron 1912 


In. 

Area 

Rad 

Long 

Tran 

Baluchistan 

IX 

511 

13 

31 - 

6 

1909 (M7) 

VII.V 

908 

17 

41 

8 


VII 

3216 

32 

60 

17 


VI 

5540 

42 

75 

24 


Table 3.4 Himalayan earthquake (1929) 


Himalayan 

1929 

(M7.1) 


So. 

Coulson 1930 

Mukherjee 1950 

In. 

Area 

Rad 

Long 

Tran 

Area 

Rad 

Long 

Tran 

VIII.V 

45973 

121 

269 

103 

85487 

165 

266 

98 

VII.V 

244421 

279 

450 

300 

339877 

329 

461 

316 

V.V 

- 


- 

- 

912236 

539 

771 

363 

III 

- 

- 

- 

- 

2226147 

842 

1134 

611 
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Table 3.5 Mach earthquake (1931) 



So. 

West 1934 


In. 

Area 

Rad 

Long 

Tran 


VIII 

2827 

30 

56 

13 

Mach 1931 

VII.V 

4536 

38 

93 

21 

(M7.4) 

VI.V 

66052 

145 

232 

75 


V.V 

113411 

190 

304 

97 


I 

301593 

310 

400 

240 


Table 3. 6 Sarigh earthquake (1931) 


Sarigh 
(1931) (M7) 

So. 

West 1 

1934 

In. 

Area 

Rad 

Long 

Tran 

VII 

531 

13 

12 

13 

V.V 

19597 

79 

111 

72 

I 

78387 

158 

231 

108 


Table 3 . 7 Quetta earthquake (1935) 


Quetta 

1935 

(M7.5) 

So. 

West] 

.936 

In. 

Area 

Rad 

Long 

Tran 

X 

1662 

23 

40 

12 

VIILV 

2733 

29.5 

52 

16 

VII 

8556 

52.2 

93 

29 

I 

256840 

286 

273 

273 


Table 3.8 Hindukush earthquake (1937) 


Hindukush 

1937 

(M7.2) 

So. 

{ 

Coulson 1938 

In. 

Area 

Rad 

Long 

Tran 

VIII 

18627 

77 

98 

60 

VII 

45239 

120 

150 . 

98 

V.V 

282743 

300 

375 

262 

IV.V 

679291 

465 

525 

413 

III 

1552603 

703 

825 

600 

II.V 

2505260 

893 

1013 

786 
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Table 3.9 Pamir earthquake (1939) 



So. 

Coulson 1940 I 


In. 

Area 

Rad 

Long 

Tran 


VIII 

5027 

40 

52 ■ 

30 

Pamir 1939 

VII 

85530 

165 

250 

108 

(M6.9) 

VI 

282743 

300 

370 

240 


V 

425447 

368 

422 

320 


IV 

622114 

445 

480 

410 


III 

817128 

510 

520 

500 


Table 3.10 Cltamba earthquake (1945) 


Chamba 

1945 

(M6.5) 

So. 

Kapila 1959 

In, 

Area 

Rad 

Long 

Tran 

X 

50 

3.98 

4 

3.8 

IX 

201 

8 

8 

8 

VIII 

9157 

54 

55 

53 

VII 

35282 

106 

106 

104 


Table 3.11 Badgom earthquake (1963) 


Badgom 

1963 

(M5.5) 

So. 

Wakhaloo 1977 

Gosain and 

Arya 1967 

In. 

Area 

Rad 

Long 

Tran 

Area 

Rad 

Long 

Tran 

VIII 

3 

0.84 

1 

0.6 

13 

2 

3 

1.5 

VII 

39 

3.48 

5 

2.5 

88 

5 

7 

4 

VI 

- 

- 

- 

- 

1357 

21 

24 

18 

V 

- 

- 

- 

- 

6685 

46 

56 

38 

IV 

- 

- 

- 

- 

22217 

84 

104 

68 


Table 3.12 Anantnag earthquake (1967) 


Anantnag 

1967 

(M5.7) 

So. 

Wakhaloo 1977 

Gosain and Arya 1 

[967 

In. 

Area 

Rad 

Long 

Tran 

Area 

Rad 

Long 

Tran 

VII 

5 

1.16 

1.16 

1 

616 

14 

22 

9 

VI 

453 

12 

12 

4 

2827 

30 

-48 

19 

V 


- 

- 

- 

4301 

37 

55 

25 

IV 

- 

- 

- 

- 

11310 

60 

80 

44 
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Table 3.13 Kinnaur earthquake (1975) 


Kinnaur 
1975 (M7) 

So. 

Hukku et al. 1975 

Singh et al. 1977 

In. 

Area 

Rad 

Long 

Tran 

Area 

Rad 

Long 

Tran 

X 

154 

7 

5 

7 

201 

8 

5 

7 

IX 

962 

17.5 

8 

19 

1018 

18 

8 

20 

VIII 

2826 

30 

20 

50 

3216 

32 

22 

50 

VII 

8821 

53 

38 

71 

9157 

54 

40 

75 


VI 

24317 

88 

72 

105 

24872 

89 

75 

108 

V 

62427 

141 

140 

135 

70650 

150 

145 

140 

IV 

160379 

226 

222 

200 

166106 

230 

232 

209 



Gosain and Arya 1967 

In. 

Area 

Rad 

Long 

Tran 

X 

2827 

30 

35 

26 

IX 

4536 

38 

46 

33 

VIII 

7854 

50 

63 

43 

VII 

11310 

60 

71 

52 

VI 

25447 

90 

105 

79 

V 

31416 

100 

122 

88 

IV 

101788 

180 

194 

172 

III 

196350 

250 

253 

243 

II 

502655 

400 

424 

388 


Table 3.14 Dharamsala earthquake (1978) 



Table 3.15 Jammu-1 earthquake (1980) 



So. 

Krishnamurthy et al. 1980 | 


In« 

Area 

Rad 

Long 

Tran 

Jammu-1 

VIII 

41 

3.6 

6 

2 

1980 (M5.5) 

VII 

114 

6 

9 

3 


VI 

314 

10 

19 

7 


V 

804 

16 

30 

9 
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Table 3.16 Jammu-2 earthquake (1980) 


Jammu-2 

1980 

(M5.4) 

So. 

Krishnamurthy etal. 1980 

In. 

Area 

Rad 

Long 

Tran 

VII 

37 

3.4 

4 

2 

VI 

255 

9 

10 

7 

V 

2641 

29 

38 

17 


Table 3.17 Dharamsala earthquake (1986) 



So. 

Gupta et al. 1986 


In. 

Area 

Rad 

Long 

Tran 

Dharamshala 

VII 

201 

8 

9 

6 

1986 (M5.7) 

VI 

1520 

22 

25 

19 


V 

4776 

39 

43 

35 


IV 

18137 

76 

126 

49 


Table 3.18 Srinagar earthquake (1988) 



So. 

Gupta 1988 

Srinagar 

In. 

Area 

Rad 

Long 

Tran 

1988 (M4.5) 

TV 

314 

10 

16 

6 


III 

2462 

28 

40 

20 


Table 3.19 Chamba earthquake (1995) 


Chamba 
1995 (M4.5) 

So. 

Pane 

le and Sharda 1995 

In. 

Area 

Rad 

Long 

Tran 

VII 

7 

1.4 

1.4 

1.4 

VI 

255 

9 

10 

7 

V 

1809 

24 

28 

21 


Table 3.20 Delhi earthquake (I960) 


Delhi 1960 
(M6) 

So. 

Mutl 

linath etal. 1969 | 

In. 

Area 

Rad 

Long 

Tran 

VII 

168 

7.3 

14 

3 

VI 

1256 ^ 

20 

22 

20 

V 

4535 

38 

49 

31 


Table 3.21 Indo-Nepal earthquake (1978) 



So. 

Kumar et al. 1981 


In. 

Area 

Rad 

Long 

Tran 

Indo-Nepal 

VI 

2123 

26 

7 

3 

1978 (M4.5) 

V 

7850 

50 

15 

7 


IV 

14520 

68 

19 

14 


III 

44466 

119 

- 

- 
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Table 3.22 Nepal earthquake (1980) 



Table 3.23 Uttarkashi earthquake (1991) 



So. 

GSI 1992 

Narula et a 

1. 1995 


In. 

Area 

Rad 

lESl 

Tran 

Area 

Rad 

Lo?5. 

Tran 

Uttarkashi 

IX 

347 

10.5 

10 

10 

565 

13 

15 

12 

1991 

VIII 

453 

12 

21 

11 

1062 

18 

26 

13 

(M6.6) 

VII 

3630 

34 

40 

33 

5730 

43 

48 

38 


VI 

26577 

92 

113 

71 

28981 

96 

123 

75 


V 

151976 

220 

265 

150 

165483 

230 

301 

175 


Narula et al. 1995 

Area 

Rad 


Tran 

415 

11 

12 

11 

939 

17 

23 

13 

5089 

40 

45 

36 

30140 

98 

123 

78 

137071 

209 

271 

161 


Table 3.24 Delhi earthquake (1994) 
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Table 3.26 Chamoli earthquake (1996) 


Chamoli 
1996 (M4.5) 

So. 

Singh and Joshi 1996 

In. 

Area 

Rad 

Long 

Tran 

VI 

452 

12 

20 

7 

V 

1257 

20 

30 

14 

IV 

5542 

42 

50 

36 


Table 3.27 Gharwal earthquake (1996) 


Gharwal 
1996 (M5) 

So, 

Pande and Ga 

rola 1996 

In, 

Area 

Rad 

Long 

Tran 

VI 

113 

6 

8 

4 

V 

452 

12 

16 

8 

IV 

3217 

32 

27 

37 


Table 3.28 Indo-Nepal earthquake (1997) 



So. 

Dasgupta et al. 200( 


Indo-Nepal 

In, 

Area 

Rad 

Long 

Tran 

1997 (M5.5) 

VI 

2827 

30 

44 

20 


V 

11310 

60 

83 

44 


Table 3.29 Chamoli earthquake (1999) 


Chamoli 

1999 

(M6.3) 

So. 

GSI 

In. 

Area 

Rad 

Long 

Tran 

VIII 

1810 

24 

40 

14 

VII 

18146 

76 

115 

50 


Table 3.30 Kathmandu earthquake (1833) 



So. 

Dasgupta 

et al. 2000 



In, 

Area 

Rad 

Long 

Tran 


IX 

16278 

72 

75 

63 

Kathmandu 

VIII 

149226 

218 

225 

200 

1833 

VI 

358727 

338 

483 

267 


VI 

520138 

407 

533 

308 


V 

970688 

556 

788 

388 


IV 

1775683 

752 

1033 

500 
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Table 3.31 Bihar-Nepal earthquake (1934) 


Bihar-Nepal 
1934 (M8.4) 

So. 

GSI 1939 

In. 

Area 

Rad 

Long 

Tran 

X 

2827 

30 

60 

15 

IX 

32669 

102 

140 

65 


70650 

150 

175 

138 

VII 

262256 

289 

256 

279 

VI 

794449 

503 

694 

368 

V 

2039858 

806 

921 

591 

IV 

3974063 

1125 

1058 

898 

II 

7302463 

1525 

1340 

1234 


Table 3.32 Bulandshaliar earthquake (1956) 



So. 

Tandon 1975 


In. 

Area 

Rad 

Long_ 

Tran 

Bulandsahar 

VIII 

1320 

20.5 

20 

17 

1956 (M6.5) 

VII 

6359 

45 

53 

38 


VI 

18618 

77 

88 

65 


V 

244421 

279 

333 

230 


Table 3.33 Roorkee earthquake (1975) 



So. 

Arya etal. 1977 

Roorkee 

In. 

Area 

Rad 

Long 

Tran 

1975 

VI 

59 

4.3 

6.5 

3 

(M4.7) 

V 

314 

10 

16 

7 


IV 

707 

15 

20 

12 


Table 3.34 Bihar-Nepal earthquake (1988) 



So. 

GSI 1993 

Thakkar et al. 1993 


In. 

Area 

Rad 

Long 

Tran 

Area 

Rad 

Long 

Tran 

Bihar- 

IX 

3019 

31 

42 

23 

- 

- 


- 

Nepal 1988 

VIII 

17671 

75 

108 

50 

35151 

106 

167 

67 

(M6.4) 

VII 

38013 

110 

160 

75 

101747 

180 

233 

139 


VI 

80425 

160 

215 

115 

189545 

246 

311 

194 


V 

113411 

190 

271 

135 

- 

- 

- 

- 
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Table 3.35 Assam earthquake (1897) 


Assam 

1897 

(M8.7) 

So. 

Narula and Sarkar 1997 

In. 

Area 

Rad 

Long 

Tran 

XI 

14527 

68 

90 

50 

X 

70686 

150 

192 

120 

IX 

180956 

240 

312 

180 

VIII 

384845 

350 

402 

300 

VII 

679291 

465 

540 

400 


So. 

Oldham 1899 

In. 

Area 

Rad 

Long 

Tran 

XI 

118177 

194 

213 

172 

VIII.V 

323549 

321 

410 

287 

VII.V 

678947 

465 

616 

418 

VI 

1719464 

740 

820 

673 

IV 

3350659 

1033 

1148 

944 


Table 3.36 Calcutta earthquake (1906) 


Calcutta 

1906 

So. 

Ml 

iddlemiss 1907 

In. 

Area 

Rad 

Long 

Tran 

VI 

12070.2 

62 

88 

44 

IV 

79382.3 

159 

190 

111 

II 

180864 

240 

278 

187 


Table 3.37 Srimangal earthquake (1918) 


Shrimangal 

1918 

(M7.6) 

So. 

Stuart 1918 

In. 

Area 

Rad 

Long 

Tran 

X 

460 

12 

15 

8 

vni.v 

4070 

36 

44 ■ 

30 

VII.V 

18377 

77 

90 

67 

VI 

40808 

114 

134 

103 

IV 

836044 

516 

744 

492 

II.V 

2371206 

869 

1020 

820 


So. 

Struat ] 

[920 

In. 

Area 

Rad 

Long 

Tran 

X 

- 

- 

- 

- 

VIII.V 

- 

- 

- 

- 

VII.V 


- 

- 

- 

VI 

22156 

84 

125 

76 

IV 

855600 

522 

750 

491 

II.V 

2215584 

840 

1018 

825 
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Table 3.38 Diihbri earthquake (1930) 


Dhubri 

1930 

(M7.1) 

So. 

Gee 1934 

In. 

Area 

Rad 

Long 

Tran 

IX 

1134 

19 

12 

15 

VII 

15829 

71 

70 

91 

VII 

62427 

141 

170 

110 

VI 

174886 

236 

273 

190 

IV 

602391 

438 

430 

381 


Table 3.39 Cachar earthquake (1984) 


Cachar 

1984 

(M5.5) 

So. 

Das 

gupta et al 2000 

In. 

Area 

Rad 

Long 

Tran 

VII 

20 

2.52 

4 

1.5 

VI 

238 

8.7 

15 

6 

V 

1018 

18 

25 

17 


Table 3.40 Assam earthquake (1950) 


Assam 

1950 

(M8.7) 

So. 


Poddar 1950 

In. 

Area 

Rad 

Long 

Tran 

XI 

6362 

45 

65 

30 

X 

90746 

170 

310 

95 

IX 

212264 

260 

460 

154 

VIII 

384650 

350 

602 

207 

VII 

638680 

451 

802 

251 

VI 

1460490 

682 

951 

351 

V 

2178812 

833 

1112 

482 

II 

4635347 

1215 

1410 

889 


Table 3.41 Assam earthquake (1975) 


Assam 

1975 

(M6.7) 

So. 

Gosavi et al. 1977 

In. 

Area 

Rad 

Lons 

Tran 

VI 

3847 

35 

59 

18 

V 

29545 

97 

131 

70 

IV 

153361 

221 

288 

174 

III 

530412 1 

411 

505 

344 
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Table 3.42 Nagaland earthquake (1970) 


Nagaland 

1970 

(M6.4) 

So. 

Guha and Gosain 1974 

In. 

Area 

Rad 

Long 

Tran 

VI 

27158 

93 

130 

70 

V 

146500 

216 

280 

155 

IV 

446286 

377 

470 

215 

in 

919019 

541 

625 

365 

II 

1650463 

725 

830 

618 


Table 3.43 Burma-India earthquake (1988) 


Burma- 
India 1988 
(M6.8) 

So. 

BIECT 

In. 

Area 

Rad 

Long 

Tran 

VII 

1256 

20 

32 

18 

VI 

6359 

45 

58 

32 

V 

20096 

85 

110 

70 


Table 3.44 Kutch earthquake (1819) 


Kutch 

1819 

So. 

Dasi 

gupta et al. 2000 

In. 

Area 

Rad 

Long 

Tran 

IX 

169007 

232 

257 

201 

VIII 

352387 

335 

389 

320 

II 

3415846 

1043 

1111 

840 


Table 3.45 Coimbatore earthquake (1900) 


Coimbat- 
-ore 1900 
(M6) 

So. 

GSI (Atl 

las) 

In. 

Area 

Rad 

Long 

Tran 

VII 

11684 

61 

80 

45 

VI 

40808 

114 

148 

84 

V 

119399 

195 

194 - 

98 

IV 

307623 

313 

274 

231 

III 

816714 

510 

510 

400 


Table 3.46 Satpura earthquake (1 938) 


Satpura 

1938 

(M6.2) 

So. 

Mukherjee 1942 

In. 

Area 

Rad 

Long 

Tran 

VII 

21632 

83 

96 

75 

VI 

120627 

196 

218 

191 

V 

278845 

298 

328 

239 

IV 

492403 

396 

464 

300 

III 

684800 

467 

546 

382 

LV 

1016610 

569 

870 

700 
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Table 3.47 Anjar earthquake (1956) 



So. 

Tandon 1956 

Dasgupta et al. 2000 



In. 

Area 

Rad 

Long 

Tran 

Area 

Rad 


Tran 


IX 

154 

7 

11 

7 

327 

10 

13 

8 


VIII 

1662 

23 

23 

24 

1806 

24 

25 

23 

Anjar 
1956 (M7) 

VII 

4070 

36 

38 

38 

4775 

39 

40 

38 

VI 

11684 

61 

72 

51 

12959 

64 

75 

55 


V 

35950 

107 

131 

90 

41186 

114 

138 

95 


IV 

110981 

188 

215 

169 

123700 

198 

225 

175 


III 

294018 

306 

343 

284 

307876 

313 

350 

280 


Table 3.48 Calcutta earthquake (1964) 



So. 

Jhingram et al. 1969 


In. 

Area 

Rad 

Long 

Tran 

Calcutta 

VIII 

11304 

60 

70 

30 

1964 

VII 

24872 

89 

95 

68 

(M5.5) 

VI 

53066 

130 

135 

89 


IV 

101736 

180 

230 

124 


II 

180864 

240 

280 

158 


Table 3.49 Koyna earthquake (1967) 


Koyna 
1967 (M7) 

So. 

Chatter ji et al. 1967 

In. 

Area 

Rad 

Long 

Tran 

VIII 

79 

5 

7 

4.5 

VII 

1134 

19 

17 

12 

VI 

1963 

25 

31 

21 


V 

133250 

206 

238 

145 

IV 

509965 

403 

463 

310 


Table 3.50 Bhadrachalam earthquake (1969) 
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Table 3.51 Broach earthquake (1970) 


'oach 

So. 

In. 

VII 

VI 

y /U 

46 1) 

V 


IV 

in 

II 



Chaubhury etal, 1970 


Area 

Rad 


- 

- 

12 

5 

1809 

24 

26 

21 

7235 

48 

60 

41 

24317 

88 

98 

79 

48281 

124 

139 

111 

102870 

181 

206 

158 


Table 3.52 Koyna earthquake (1973) 


Koyna 

1973 

(M5.2) 


Table 3.53 Shimoga earthquake (1975) 



So. 

Gosavi et a 

11977 

In. 

Area 

Rad 

Lonj 

V 

616 

14 

16 

IV 

3317 

32.5 

38 

III 

18137 

76 

91 

II 

88624 

168 

205 





So. 


In. 

Shimoga 

V 

1975 (M5) 

IV 


III 


II 



Gosavi etal. 1977 


Area Rad Lon 


7235 


28939 


76416 


Table 3.54 Jaisalmer earthquake (1991) 


Jaisalmer 

1991 

(M5.5) 


Table 3.55 Saurashtra earthquake (1993) 



So. 

Dasgupta et al. 2000 

In. 

Area 

Rad 

Long 

Tran 

VIII 

114 

6 

3 

3 

VII 

908 

17 

45 

6 

VI 

13267 

65 

96 

37 

V 

36625 

108 

129 

88 
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Table 3.56 Killari earthquake (1993) 


Killari 

1993 

(M6.2) 


So. 

GSI (Atlas) 

In. 

Area 

Rad 

Long 

Tran 

IX 

154 

7 

5 

3 

VIII 

453 

12 

9 

5 

VII 

616 

14 

16 

11 

VI 

8168 

51 

57 

38 

VI 

14520 

68 

148 

100 


Table 3.57 Bonaigarh earthquake (1995) 


Bonaigarh 
1995 (M4.4) 

So. 

Dasj 

gupta et al. 2000 

In. 

Area 

Rad 

Long 

Tran 

V 

1520 

22 

19 

23 

IV 

6645 

46 

47 

41 

III 

15829 

71 

76 

65 


Table 3.58 Jabalpur earthquake (1997) 



So. 

Rai et al. 1997 


Dasgupta et al. 2000 


In. 

Area 

Rad 

Long 

Tran 

Area 

Rad 

Long 

Tran 

(abalpur 

VIII 

320 

10 

12 

9 

360 

12 

14 

11 

997 (M6) 

VII 

2262 

27 

36 

20 

2642 

29 

42 

18 


VI 

8708 

53 

84 

33 

10568 

58 

92 

35 


V 

74374 

154 

178 

133 

84496 

164 

188 

143 


Table 3.59 Bhuj earthquake (2001) 


Bhuj 2001 
(M7.9) 

So. 

Narula and Chaubey 2001 

GSI 2( 

)01 

In. 

Area 

Rad 

Long 

Tran 

Area 

Rad 

Lon^ 

Tran 

X 

2324 

27 

42.5 

18 

1005 

18 

20 

16 

IX 

11011 

59 

90 

42 

11517 

- 61 

78 

47 

VIII 

30049 

98 

170 

85 

41639 

115 

141 

94 

VII 


- 

- 

- 

104458 

182 

266 

125 
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Table 3. 60 Earthquake details of Northwest province 


SI. No 

Name 

Date 

Mag, 

Foe. Dep. 

Max. int. 

Co-ordinate 

Source 

B 

Kashmir 

30-May-1885 

7 

12 

1/MS 

34.1 

74.61 

Jones 1885 

VlII/MSK 

Gosain and Aiya 1 967 

2 

Kangra 

4-Apr-l905 


32.6-62 

X/MMl 

33.6 

76.2 

Middlemiss 1910 

3 

Baluchistan 

21 -Oct- 1909 

- 

- 

IX/MMI 

30 

68 

Haron 1910 

4 

Himalaya 

1 -Feb- 1929 

7 1 

200-250 

VIII/RF 

36 5 

70.5 

CoLiIson 1930 

VIII/RF 

Mukei]ee 1950 

5 

Mach 

25-Aug-1931 

7.4 

- 

VIII/RF 

30.2 

67.7 

West 1934 

6 

Sharigh 

27-Aiig-1931 

7 

- 




West 1934 

7 

Quetta 

31-May-1935 

7.5 

- 



66.8 

West 1936 


Hindukush 

14-NOV-1937 

WHi 


VIII/RF 

37.5 

71 5 

Coulson 1938 

BB 

Pamir 

21-NOV-1939 





— 

74 

CoLilson 1940 

10 

Chamba 

22-Jul-1945 


60 




Rapila 1959 

11 

Badgom 

2-Sep-1963 

5.5 





Wakhaloo 1977 

12 

Anantnag 

20-Feb-1967 

5.7 

24 

VI I/M MI 

33.7 



75.28 

Wakhaloo 1977 

5.5 

20 

VII/MMI 

Gosain and Arya 1967 

13 

Kinnaur 

19-Jan-1975 



X/MMI 

32.5 

— 

78.43 

Gosavi et al 1 977 


X/MMI 

Singh et al 1 977 

7 


X/MMI 

Hukku<?/^// 1977 

14 

Dharamsala 

1 5-Jun-1978 

5 

- 


32.2 

76.61 

Kumar et al 1 98 1 

15a 

Jammu-1 

24-ALig-1980 

5.5 

- 

VIII/MMI 

- 

- 

Knshnamurthy er al 

1990 

16 

Jammu-2 

24-Aiig-1980 

5.4 

- 

VII/MMI 



Knshnamurthy ei 7/ 

1990 


Dharamsala 

26-Apr-I986 


23 

VIII/MMI 

32.3 

76.3 

Gupta 1986 


Srinagar 

8-Feb-1988 

3.2 

- 

IV/MSK 

34.4 

75.3 

Gupta 1988 


Chamba 

24- Mar- 1995 

■Bl 

10 

VI/MSK 

- 

- 

Pande and Sharda 1995 


Table 3.61 Earthquake details of Northern province 


SI. No 

Name 

Date 

Magnit 

tide 

Focal 

depth 

Max. inten 

Co-ordinate 

Source 

20 

Delhi 

27-Aug-1960 

6 

- 

VII/MMI 

28.2 

77.4 

Muktinath et al 1 969 

21 

Indo-Nepal 

21 -May- 1979 

5.8 

13-16 

IV/MMI 

30 

80.3 

Kumar err//. 1981 

22 

Nepal 

29-Jul-1980 

6.1 

- 

VIIl/MMl 

29.6 

81.07 

Snvastava et al 1980 

23 

Utterkashi 

20-Oct-I991 

6.6 


VIll/MMI 



Thakkar 1993 

- 

IX/MSK 

- 

- 

GSI 1992 

^6.4 

- 

IX/MSK 

- 

- 

Narulae/r/Z 1995 

24 

Calcutta 

15-Apr-1964 

5.5 

36 

VIII/MMI 

21 

88.5 

Jhingran et al 1969 


Delhi 

28-Jul-1994 ^ 

4 

4 

V/MSK 

- 

- 

Gupta and Sharda 1 996 


Bhind 

1 -Sep- 1994 

4.8 

- 

7/MSK 

- 

- 

Pande err// 1995 

27 

Chamboli 

23-Jan-1996 

4.5 

20+-5 

VI/MSK 

- 

- 

Singh and Joshi 1996 

28 

Garhwal 

26-Mar-l996 

5 

- 

Vl/MSK 

- 

- 

Pande and Gairola 

1996 

29 

Indo-Nepal 

5-Jan-1997 

5.5 

16 

VI/MSK 

29.8 

80.5 

Kumar err//. 1981 

30 

Chamboli 

29-Mar- 1999 

6.3 

21 

Vlll/MSK 

30.5 

79.4 

Dasgupta et al. 2000 
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Table 3.62 Earthquake details of Foredeep province 


31 

Kathmandu 

26-Aiig-l833 

- 

- 

IX/MMI 

- 

- 

DasgLipta et a! 2000 

32 

Bihar-Nepal 

]5-Jan-1934 

8.4 

20-30 

X/MMI 

26.5 

86.5 

Gsi 1939 

33 

Bulandshahar 

lO-Oct-1956 

6.5 

20 

VIII/MMI 

28.2 

77.8 

Tandon 1975 

34 

Roorkee 

6-NOV-1975 

4.7 

- 

VI/MMI 

32.4 

78.5 

Arya et al 1977 

35 

Bihar-Nepal 

21 -Aug- 1988 

64 

65 

IX/MMl 

26.7 

86.63 

GSI 1993 


Table 3.63 Earthquake details of Northeast province 


SI. No 

Name 

Date 

Mag. 

F.D. 

Max. int. 

Co-ordinate 

Source 

36 

Assam 

12-JLnvl897 

8.7 

155-310 

XII/MMI 

26 

91 

Oldham 1899 

Narula and Sarkai 

1997 

37 

Calcutta 

29-Sep.-06 

- 

- 

VI/MMl 

27 

97 

Middlemiss 1907 

38 

Shrimangal 

S-Jiil-1918 

7.6 1 

12.4-14 

X/MMI 

24 5 

91 

Staint 1918 

39 

Dhubri 

3-Jul-1930 

7.1 

. 

IX/MMI 

25.5 

90.9 

Gee 1934 

40 

Cachar 

31 -Dec- 1984 

5.5 

2 

VII/MMI 

24.7 

92.85 

Dasgupta et al 2000 


Table 3.64 Earthquake details of Subduction province 


41 

— 

Assam 

15-Aug-1950 

8.7 

46.5 

X/MMI 

28.5 

96.5 

Poddar 1950 

X/MMI 

Dasgupta et al 2000 

42 

Assam 

S-Jul-1975 

67 

60 

VII/MMI 

21.5 

94.7 

GQSzy\etaI 1977 

43 

Nagaland 

29-Jul-1970 

6.4 

59 

VII/MMI 

26 

95.4 

Guha and Gosavi 

1974 

44 

Burma-India 

6-Aug-1988 

6.8 

- 

VII/MMI 

. 

- 

BIEST 
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Table 3. 65 Earthquake details of Peninsular province 


45 

Kutch 


- 

- 

IX/MMl 

- 

- 

Dasgupta et cil 2000 

46 

Coimbatore 

8-Feb-1900 

6 

70 

VII/MMI 

10.8 

76.75 

Basu 1 964 

47 

Satpura 

14- Mar- 1938 

6.2 

43 

Vll/MMI 

21.53 

75.83 

Mukhei]ee 1942 

48 

Anjar 

21-Jul-1956 

7 

13-18 

IX/MMI 

23.3 

70 

Tan don 1959 

IX/MMI 

Dasgupta et al 2000 

49 

Koyna 

1 1- Dec- 1967 

5 

1 

VIII+/MM1 

17.2 

73.72 

Chatterji 1967 

50 

Bhadrachalam 

13-Apr-1969 

6.5 

40 

Vil/MMI 

17.9 

80.6 

Mukherjee 1971 

5.3 

25 

VII/MMI 

17.8 

80.67 

Dasgupta et al 2000 

51 

Broach 

23-Mar-70 

6.7 

10 

VII/MMl 

21.68 

73.2 

Chaudhury et al 1970 

8 

VII/MMI 

Dasgupta et al 2000 

52 

Koyna 

17-Oct-1973 

5.2 

10 

VI/MMI 

17.7 

73.9 

Gosavierr// 1977 

53 

Shimoga 

12-May-1975 

5 

35 

V/MMI 

13.8 

75.3 

Gosnvietal 1977 

54 

Jaisalmer 

8-Sep-1991 

5.5 

30 

VlII/MMI 

26.4 

70.66 

Dasgupta et al 2000 

55 

Saurashtra 

24-ALig-1993 

- 

29 

V/MMI 

20.7 

71.44 

Dasgupta et al 2000 

56 

Killari 

30-Sep-1993 

6.2 

10 

IX/MMI 

18.1 

76.5 

Dasgupta et al 2000 
Narulaera/ 1996a 
Mohan and Rao 1 994 
Ramaknshna et al 

1996 

57 

Bonaigarh 

27-Mar-1995 

4.4 

21 

V/MSK 

21.7 

85.59 

Dasgupta et al 2000 

58 

Jabalpur 

30-Mar- 1900 

6 

36 

Vm/MSK 

23.1 

80.12 

Dasgupta (T// 2000 
ChatLirvedi et al 2000 

59 

Bhuj 

26-Jan-2001 

7.9 

- 

X/MSK 

23.8 

70.3 

Nani la 

GSI 2001 
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Figure 3.3 Isoseismal of Baluchistan earthquake of 21^' October 1906 following RF scale 
after Haron (1910). 
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Figure 3.4 Isoseismal of Northwest Himalayan earthquake (M 7.1) of f ‘ February 1929 
following RP scale after Coulson (1930). 
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Figure 3.5 Isoseismal of Northwest Himalayan earthquake (M 7. 1 ) of 1 Februarv 1 929 
following MM scale after Mukherjee (1950). 
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Figure 3.6 Isoseismal of Mach earthquake (M 7.5) of 27'’^ August 1931 following RF scale 
after West (1934). 
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Figure 3.7 Isoseismal of Sharigh earthquake (M 7) of 25'" August 1931 following RF scale 
after West (1934). 




Chapter 3 Review of Ecirthqiitikc Intensity Database 


100 


- 7/ y/ 

'P 














x.vckMj. 













OAISALMER 


^:^r!^'V''-‘^s , :fk-r''-*-\\ •y^.-'.i^'^^ )\’'^lL ^ ‘X* 
yat l ■■•'•.■VvVr '^-''. w.kaA^»“'4>--~-^:^^ 






^‘<0. >, X^\ f jViie^Mu N>w‘^-*«') \ 


Figure 3.7 Isoseismal of Sharigh earthquake (M 7) of 25'*’ August 1931 following RF scale 
after West (1934). 
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Figure 3.8 Isoseisnial of Quetta earthquake (M 7.5) of 31^* May 1935 following RF scale 
after West (1936). 
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Figure 3.9 Isoseismal of Hindkush earthquake (M 7.2) of 14'*^ November 1937 following RF 
scale after Coulson (1938). 
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Figure 3.10 Isoseismal of Pamir earthquake (M 6.9) of 21®' November 1939 following MM 


scale after Coulson (1940). 
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Figure 3.13 Isoseismal of Anantnag earthquake (M 5.7) of 20'*’ February 1967 following 
MM scale after Wakhaloo (1977). 
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Figure 3.14 Isoseismal of Anantnag earthquake (M 5.7) of 20*'^ February 1967 following 
MM scale after Gosain and Arya (1967). 
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Figure 3.15 Isoseismal of Jammu earthquake (M 5.5, 5.4) of 24*'' August 1980 following 
MM scale after Krishnamurthy et al (1990). 
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Figure 3.16 Isoseismal map of Srinagar earthquake (M 3.2) of 8'*’ February 1988 following 
MSK scale after Gupta (1988). 
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Figure 3.17 Isoseismal map of Chamba earthquake (M 4.5) of 24"' March 1995 following 
MSK scale after Pande and Sharda (1995). 
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Figure 3.18 Isoseismal map of Kathmandu earthquake of 26''^ August 1833 following MM 
scale taken from Dasgupta et al. (2000). 
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Figure 3.21 Isoseismal map of Bulandshahar earthquake (M 6.5) of 10''’ October 195< 
following MM scale after Tandon (1975). 
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Figure 3.22 Isoseismal map of Delhi earthquake (M 6) of 27'’’ August 1956 following MM 
scale after Muktinath et al. (1969). 
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Figure 3.24 Isoseismal of Kinnaur earthquake (M 7) of 19''' January 1975 following MM 
scale after Gosavi et al. (1977). 
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Figure 3.25 Isoseismal of Kinnaur earthquake (M 7) of 19*'’ January 1975 following MM 
scale after Singh et al. (1975). 
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Figure 3.27 Isoseismal map of Roorkee earthquake (M 4.7) of 6'*' November 1975 
following MM scale after Arya et al. (1977). 
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Figure 3.29 Isoseismal map of Nepal earthquake of (M 6.1) 29^'’ July 1980 following MM 
scale after Srivastava et al. (1980). 
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Figure 3.31 Isoseismal map ofBihar-Nepal earthquake (M 6.4) of 20*’' August 1988 
following MM scale after GSI (1993). 
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Figure 3.32 Isoseismal map of Bihar-Nepal earthquake (M 6.4) of 20*'^ August 1988 
following MM scale after Thakkar (1988). 
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Figure 3.33 Isoseismal map of Uttarkashi earthquake (M 6.6) of 20'*^ October following 
MSK scale after Narula et al. (1995). 



Chapter 3 Review of Earthquake Intensity Database 


127 



Figure 3.34 Isoseismal map ofUttarkashi earthquake (M 6.6) of 20'*’ October following 
MSK scale after Rastogi and Chadha (1995). 
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Figure 3.37 Isoseismal map of Bhind earthquake (M 4.8) of September 1994 followin 
MSK scale after Pande et al. (1995). 
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Figure 3.39 Isoseismal map of Garhwal earthquake (M 5) of 26‘'^ March 1996 following 
MSK scale after Pande and Gairola (1996). 





C 'hapter 3 Review oj Ecu thquake Intensity Database 


13 



Figure 3.40 Isoseismal map of Indo-Nepal earthquake (M 5.5) of 5*'^ January 1997 followin 
MM scale after Dasgupta et al. (2000). 
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Figure 3.42 Isoseismal map of Assam earthquake (M 8.7) of 12'’’ June 1897 following 
Olaham's scale after Oldham (1899). 
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Figure 3.43 Isoseismal map of Assam earthquake (M 8.7) of 12*'’ June 1 897 following MSK 
after Narula and Sharda (1997). 
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Figure 3.44 Isoseismal map of Calcutta earthquake of 29'*’ September 1906 following MM 
scale after Middlemiss (1907). 
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Figure 3.45 Isoseismal map of Srimangal t 
MM scale after Staurt (1918). 
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lake (M 7.6) of 8* July 191 8 following 
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Figure 3.47 Isoseismal map of Assam earthquake (M 8.7) of 15"’ August 1950 following 
MM scale after Poddar (1950). 
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MM scale taken from Dasgupta et al. (2000). 
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Figure 3.51 Isoseismal map of Cachar earthquake (M 5.5) of 31®' December 1984 followin 
MM scale after Narula [taken from Dasgupta et al. (2000)]. 
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Figure 3.52 Isoseismal map of Buima earthquake (M 6.8) of 6"' August 1988 following MM 
scale after Kumar (1992). 
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Figure 3.53 Isoseismal map of Kutch earthquake of 16’*’ June 1819 following MM scale 
taken from Dasgupta et al. (2000). 
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Figure 3.54 Isoseismal map of Coimbatore earthquake (M 6) of 8th February 1900 

following MM scale after Basil (1964) taken from Kaila and Sarkar (1978). 
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Figure 3.55 Isoseismal map of Satpura earthquake (M 6.2) of 14*'’ March 1938 following 
MM scale after Mukherjee (1942) taken from Dasgupta et al. (2000). 
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Figure 3.56 Isoseismal map of Anjar earthquake (M 7) of 21*‘ July 1956 following MM 
scale after Tandon (1959). 






Figure 3.57 Isoseismal map of Anjar earthquake (M 7) of 21^‘ July 1956 following MM 
scale taken from Dasgupta et al. (2000). 
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Figure 3.58 Isoseismal map of Calcutta earthquake (M 5.5) of 15'’’ August 1964 following 
MM scale after Jhingran et al. (1969) taken from Kaila and Sarkar (1978). 
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Figure 3.60 Isoseismal map of Bhadrachalam earthquake (M 6.5) of 13"^ April 1969 
following MM scale after Mukheijee (1971) taken from Kaila and Sarkar 
(1978). 







Figure 3.62 Isoseismal map of Broach earthquake (M 6.7) of 23'^’^ March 1970 following 
MM scale after Chaudhury et al. (1970) taken from Kaila and Sarkar (1978). 
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Figure 3.64 Isoseismal map of Koyna earthquake 
MM scale after Gosavi et al. (1977). 


5.2) of 17''^ October 1973 following 
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Figure 3.65 Isoseismal map of Shimoga earthquake (M 5) of 12“’ May 1975 following MM 
scale after Gosavi et al. (1977) taken from Kaila and Sarkar (1978). 
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Figure 3.66 Isoseismal map of Jaisalmer earthquake (M 5.5) of 8^'’ November 1991 
following MM scale taken from Dasgupta et al. (2000). 
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Figure 3.67 Isoseismal map of Saurastra Coast earthquake of 24"’ August 1993 followin, 
MM scale taken from Dasgupta et al. (2000). 
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Figure 3.68 Isoseismal map of Killari (Latiir) earthquake (M 6.2) of 30'" September 1993 
following MM scale taken from Dasgupta et al (2000). 
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Figure 3.69 Isoseismal map of Killari (Latur) earthquake (M 6.2) of 30*'’ September 1993 
following MM scale after Narula et al. (1996a). 
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Figure 3.70 Isoseismal map of Bonaigarh earthquake (M 4.4) of 27''^ March 1995 following 
MSK scale taken from Dasgupta et al. (2000). 
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Figure 3.71 Isoseismal map of Jabalpur earthquake (M 6) of May 1997 following MSK 
scale taken from Dasgupta et al. (2000). 
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Figure 3.72 Isoseismal map of Jabalpur earthquake (M 6) of May 1997 following MSK 
scale after Chaturvedi et al. (2000). 
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Figure 3.73 Isoseismal map of Bhuj earthquake (M 7.9) of 26'*' January 2001 following 
MSK scale after Narula and Chaubey (2001). 
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Figure 3.75 Plot of earthquakes considered in the present study 




Chapter 4 

Development of Attenuation Relationships 

4.1 Introduction 

This chapter studies the attenuation of intensity with distance from the intensity- 
distance data obtained from past 59 earthquakes from all over India. Depending on the 
regression analysis results and general geology of the region, Indian subcontinent is 
divided into six attenuation provinces namely Northeast, North, Northwest, Peninsular, 
Foredeep and Subduction province. New attenuation relationships are developed for each 
province and they are compared with the existing attenuation relationships available in 
the literature. In addition, intensity attenuation relations are developed for whole Indian 
subcontinent. 

4.2 Modeled Shape of Isoseismal 

In the literature most of the equations are developed assuming that the isoseismals 
represent a family of concentric circles. Radius of a particular isoseismal is taken as 
epicentral distance for coiresponding intensity isoseist. But the shape of isoseists is 
elongated along the length of the fault, which gives rise to considerable difference 
between the attenuation along the major and minor axis of the isoseismals (Chandra et ciL. 
1979). Considering this point, two types of models are used in the present work. One 
representing a family of concentric circles and the other representing a family of 
concentric ellipses. In the first case, attenuation relations are developed for the radius of 
the each intensity circles while in the second case intensity attenuation lelationships aie 
developed along the major and minor axis. 

4.3 Attenuation Relationships of Epicentral Intensity 

form of equation depends whether intensity is proportional to the logarithm of 

seismic energy density, £, at any location or is proportional to a power of the energy 
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density. Both approaches are used in the present study. The type of attenuation relation 

depends on how the intensity is related to the seismic energy density. If it is assumed that 
the intensity is proportional to the logarithm of seismic energy density {£) at any location, 
the attenuation relationships can be written as, 

1 = !<! +kjn E. (4.1) 


where kj and Ay, are the proportionality constants. 


Using the classical equation of energy decay 
A £ 

^ n 


£ = 


^0 


4n 


(4.2) 


where E,, is the total energy released, ‘Vf’ is a constant representing the geometric 

spreading, and “c” is a constant representing rate of absorption, and A is the hypocentral 
distance. From Equations 4.1 and 4.2, 


E \ 

I = k_,+k,ln —^\-k^d{lnA)-k^cA 
[4nJ 


(4.3) 


where k^ , k^ are the proportionality constants and A is the focal depth. 


Using the abbreviations 

a I = AjC {inA) + k ,(l A 
d I - k^d 
c, =k,c 


and eliminating E„ using Equation (4.2), we get 

I = I ^ + a, -c,A- d ilnA (4.4) 

If the effect of exponential absorption is small compared to that of geometric 
spreading, this becomes equivalent to equations used by Gutenberg and Ricther (1942) of 
the fonn 


/ = /„ +< 7 y -diltiJ . 


(4.5) 
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H geometric spreading is neglected as compared to absoiption, the equation reduces to 
l=l„+a,-CiA (4.6) 

One of the simplest possible alternatives to Equation 4.4 is that intensity 
proportional to a power of the energy density 

I = kE^ ( 4 . 7 ) 

which combined with Equation 4.2 gives 


I = k 


(EA'‘ 


'0 

\4TtJ 


(Ip -(/) I 




At the epicenter 


=k 


y4Ttj 






Dividing Equation 4.8 by Equation 4.9 
I = I 


(4.8) 


(4.9) 


(4.10) 


Taking the logarithm of both sides of this and using the abbreviations 

a, = cl p InH + c p H 

cL = cl p 

c.-cp 


we get 

Inl = Inl ^ + a. - c, A - dpiA (4.11) 

If di is very small, it becomes 

Inl = lnp+ci,-C:A (4.12) 

If Cl: is also very small 

lnI = lnp-d,A (4.13) 

If a 2 and C} are both very small 


In I = Inl - cl Jn A 


(4.14) 



Chapter 4 Development of A twmiation Relationships 


172 


The difference between Equations 4.4 and 4.11 lies in the curvature of the 
variation of intensity with distance. Taking derivative of Equation 4.4 
dl dj 


dA 


zJ 


(4.15) 


Taking the derivative of Equation 4,1 1 


dJ ~ 


cl. 

— + c 
^ J 


(4.16) 


Since intensity 7’ is itself a decreasing function of A, the product of I and 


zl 


■ + c. 


shows that Equation 4.11 is steeper near the epicenter and flatter at larger 


distances than Equation 4.4. 

All of the equations describe families of parallel curves whose height above the 
abscissa depends only on/„. Inspection of isoseismal maps suggests that attenuation of 

intensity with distance is often faster for small earthquakes, so that the family of lines 
should diverge. To consider earthquake size effect on attenuation, two equations are 
proposed. 


lnI=lnI„rC,{A~A,„)-dJn — 


(4.17) 


The intensity attenuation equations are based on the assumption that the energy is 
radiated from the point source through a space of simple geometry (e.g., a unifomi 
hemisphere only a function of depth). The equations do not give good results when 
distances are not large when compared to the source dimensions. An alternative approach 
would be to postulate that the intensity-distance relation is the sum of two attenuation 
rates, one more important at short distances, the other at large distances. A near-source, 
rapidly decreasing term should be added to the trend obseiwed at large distances. This will 
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take the form of the attenuation at large distance being related to energy (magnitude) of 

the earthquake, and that the near-source being more closely linked to maximum intensity, 
that includes allowance for many factors besides energy. A simple relation of this sort is 

7 = //?Z(-CjZ) + dj— (4.18) 

■ ^ A 

where M is the magnitude. In this equation cijM has been substituted for /„ and chh/A 
replaces the constant cij. The r/j term provides for a departure from the general form of 
Equation 4.3 at short distances. 

Another family of equations similar to Equations 4.4 to 4.10 can be obtained by 
substituting {1.5M-1.5) for /„ based on the Gutenberg and Richter’s (1956) empirical 

relation between magnitude and MM intensity 

(''•19) 

3 

The equation will look like 

I = a->rb M + c A + d log, ^{ a) 

It can be assumed that magnitude is directly related to the maximum intensity 
(Lee and Trifunac, 1985), hence new attenuation relationships can be derived involving 
the maximum intensity as the earthquake parameter 

I ^ Cl +b +cA + d Log,o (d) (4.20) 

In all of the above equations, the coefficient “tf ’ defines the rate of geometric 
spreading of the energy, the coefficient “c” defines the rate of exponential absorption and 
and “/i” tenns are related to the boundary conditions at the source. Instead of the 
hypocentral distance, epicentral distance can also be used due to the lack of focal depth 
data. In the present study, both distances are used and for events whose focal depth are 
not available, same are suitably assumed to be equal to the average focal depth of the 



Chapter 4 Development of Aitemiation Relationships 


174 


region concerned. The equations involving the logarithm of intensity are not reported in 
the present study because they are found to be not good for Indian subcontinent [Chandra, 
1982], 

4.4 Attenuation Provinces 

Following the general geology of the region and the standard deviations obtained 
from the preliminary analysis (Section 4.5.8) the Indian subcontinent is divided into six 
attenuation provinces viz. Northeast, Northern, Northwest Peninsular, Foredeep and 
Subduction provinces which conforms to the attenuation provinces of Kaila and Sarkar 
(1978) and Chandra (1992). Attenuation relationships have been separately analyzed and 
by combining all the provinces, a single equation for the whole Indian subcontinent is 
also derived. 

4.5 Development of Attenuation Relationships 

4.5.1 Selection of Parameters 

The parameter, which is being predicted, is known as dependent variable and the 
parameters used for the predictions are called independent variables. In the present study, 
intensity on MM scale is taken as the dependent variable. 

4.5.2 Earthquake Parameters 

Seismic magnitude is one of the ideal parameters, which can be used to lepiesent 
the earthquake size. But for past earthquakes, the magnitudes estimated are very 
approximate and are liable to be incorrect because of lack of instrumental data. Keeping 
this in mind, equations are developed considering other parameters like the epicentral 
intensity, maximum intensity and a combination of those. Maximum intensity is not 
frequently used as attenuation parameter (Lee and Trifunac, 1985) but it is good foi 
practical purposes since it directly gives the damage. In literature most of the equations 
are based on the epicentral intensity, but is a fictitious quantity depending upon the 
method employed to estimate it and will differ from method to method. Hence for a 
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particulai earthquake, what would, be the exact epicentral intensity is a question and only 

the drop of intensity with distance can be estimated using the derived equation with 
epicentral intensity as the earthquake parameter. Due to the above reason the equation 
involving the epicentral intensity is not good in predicting the intensity for a future 
earthquake, in spite of smaller standard error. This equation gives the drop of intensity 
with distance and the drop is only due to distance parameter, while the other models have 
dual dependency. Hence it gives true values of the absoiption and geometric spreading 
constants, independent of the earthquake size. In the present study, different earthquake 
parameters like magnitude, maximum intensity and the epicentral intensity have been 
used. In addition to that, equations are developed to estimate the epicentral intensity for a 
given magnitude. 

4.5.3 Independent Variables 

In the present study, equations are developed by considering different variables 
like magnitude of earthquake, M, maximum mapped intensity, /„,«\ or epicentral intensity 
/„. Some equations are developed by taking more than one independent variable like 
magnitude and maximum intensity simultaneously. 

4.5.4 Source Mechanism 

In the present study, isoseismals from 59 earthquakes are used as data set. Since 
fault mechanism of all the events are not well defined, hence the style of faulting is not 
included as a variable in the present study. 

4.5.5 Propagation Parameters 

The propagation parameters characterize the effects due to wave scattering, 
inelastic attenuation and geometric attenuation of intensity. Since the type and length of 
faulting, amount of slip etc. are not available, the distance measurements such as distance 
to horizontal projection of fault rupture, closest distance to rupture, distance to zone of 
maximum energy released etc. could not be used as the propagation parametei. In this 
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study epicentral distance as well as the hypocentral distance are used as the propagation 
parameters. 

4.5.6 Site Parameters 

Other than the broad classification of the site conditions [Lee and Trifunac, 1985], 
no detailed information is available about the local site geology. The intensity isoseists 
are very rough and due to the lack of information and approximation the site parameters 
are not included in the present study. 

4.5.7 Type of Equation Used 

There are many types of attenuation relationships existing in the literature 
involving epicentral intensity and magnitude (Howell and Schultz, 1975). Different types 
of equations involving magnitude and epicentral intensity have been used in the present 
study. Some new types of intensity attenuation relationships involving the maximum 
intensity and combinations of maximum intensity and magnitude have been tried, which 
have resulted in good correlation. The types of equation used are 
Equation 1 
I = a + b, M + c R 
Equation 2 

I = a + b, M + d log I,) {r) 

Equation 3 

I = a + b, M + c R + d log^oi^) 

Equation 4 
l = a + b, + cR 


Equation 5 

I = a + b, -fc//og;o(i?) 
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Equation 6 

I = a+b, + cR + dlogM 

Equation 7 

I = Cl + b, +h.M + cR + d log,, (r) 

Equation 8 

/-/„ =a + cR 

Equation 9 

/-/„ =a + d log{R) 

Equation 10 

I -I ~a + c R + d log{R). 

In the above equation the constant “c” defines the rate of absorption and the 
constant “rf ’ defines the loss due to geometric spreading. 

4.5.8 Preliminary Analysis of the Data 

Before the data is regressed to get a particular attenuation model, it is necessary to 
do some preliminary analysis to eliminate inconsistent data as well as the foiTn of 
variables that gives a better fit for the data. The regression analysis is carried with the 
above models and standard deviation all each data points are calculated. The data points 
with relatively larger standard deviation are removed from the data set and not considered 
for the final regression analysis. 

4.6 Development of Intensity Attenuation Relationships 

Tables 4.1 to 4.20 give the coefficients of regression related to each region along 

with their confidence band, t-statistic values, cr and percentage variation explained values 
{R-) of the model. Separate regression analysis has been done for the whole India, to have 
a single equation to describe the attenuation of whole Indian subcontinent. Earthquakes 
are grouped in two categories one with magnitude equal to or more than 6 and the other 
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with magnitude less than 6. The two categories have been separately analyzed keeping in 
mind that the smaller earthquake attenuates at much faster rate than that of larger 
earthquake hence, they have to be separated out. 

4.7 Adequacy of the model 

The adequacy of the model can be checked either by analyzing the residuals or by 
testing the significance of each coefficient derived during the analysis. 

4.7.1 Analysis of Residue 

The overall adequacy of the model is best assessed from an analysis of residuals 
By the veiy nature of regression analysis, the residuals have a mean close to zero. If no 
apparent trend in the residuals is observed, then the model can be considered adequate. A 
trend indicates an inadequacy in the model in the model to predict the data and would 
require modification of the functional form. The residuals are plotted with respect to 
epicentral distance {R) as shown in Figures 4.7 to 4.10. No particular trend is shown in 
any plots, from which it can be ascertained that the models used in the regression analysis 
are adequate to represent the variation in data. 

4.7.2 Significance Testing of Coefficients 

The significance of the calculated coefficients can be estimated by means of a t- 

test. The r-statistic is defined as 

t = 

SE[i,) 

Under the nonnality assumption, the variable b, follows t distribution with N-P degree of 
freedom, where P is the number of coefficients in the equation. If the true value of hi is 
specified in the null hypothesis, the t-statistic can be computed from the available sample 
and used as a test statistic. 

All the coefficients are tested with a level of significance of 5% for the particular 
value being equal to zero. If the null hypothesis of these coefficients being equal to zeio 
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can not be rejected when tested with a level of significance of 5%, then the values for 
these events should be used with caution. 

4.8 Squared Multiple R 


The squared multiple R represents the proportion of variation in the dependent 

variable, accounted for by the combination of independent variable. It is also called as the 

multiple correlation coefficient and defined as 

_ {Sittit of Squares due to regression) 

{Total sum of squares, corrected for the mean Y ) 




where i is the number of data points, is the predicted dependent variable, Y is the mean 


and Y is actual value obseiwed. The value of R' lies in between 0 to 1 . A perfect fit of the 
data between f and Y would give R' =1, while values closer to one suggest a good fit of 


the data for a given model. For the present study the R' values are calculated for each 
case and used to detennine whether the particular model is good or not. 
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Table 4.9 Values of constants, confidence band, t-statistic values, standard deviation, standard error of the equation and values of the model (equation 3 

with hypocentral distance). 
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Table 4.10 Values of constants, confidence band, t-statistic values, standard deviation, standard error of the equation and values of the model (equation 

3 with epicentral distance). 
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Table 4.17 Values of the regression constants of the equation Ij,,^^=e 4 -fM 
considering all earthquakes of magnitude more than or equal to 6, 


Province 

e 

Ge 

f 


G 

Northern 

1.328 

1.38 

1.070 

0216 

1.100 

Northwest 

1.729 

0.809 

0.947 

0.138 

0.804 

Northeast 

-4.011 

2.294 

1.718 

0.336 

1.050 

Peninsular 

-0.670 

1.730 

1.345 

0.291 

1.010 

Indian sub. 

0.526 

0.839 

1.125 

0.133 

1.006 


Table 4.18 Values of the regression constants of the equation I,„^^=e + fM 
considering all earthquakes of magnitude less than 6. 


Province 

a 

Ge 

b 

CJf 

G 

Northern 

1.711 

1.807 

0.913 

0.350 

0.706 

Northwest 

0.854 

2.795 

1.158 

0.530 

0.874 

Peninsular 

-4.733 

7.425 

2.176 

1.501 

0.832 

Indian sub. 

0.359 

1.662 

1.198 

0.323 

0.7198 


Table 4.19 Values of the regression constants of the equation 1^ =e + f M + 


considering all earthquakes of magnitude more than or equal to 6. 


Province 

e 

Ge 

f 

CTf 

g 



a 

Northern 

-4.430 

1.399 

0.321 

0.170 

1.403 

0.196 

0.83 

Northwest 

5.776 

2.217 

1.41 

0.338 

-0.275 

0.175 

1.22 

Northeast 

-2.728 

2.364 

1.203 

0.610 

0.436 

0.315 

0.88 

Peninsular 

2.183 

2.557 

0.531 

0.520 

0.450 

0.279 

1.37 

Indian sub. 

1.472 

1.000 

0.202 

0.270 

0.704 

0.140 

070 


Table 4.20 Values of the regression constants of the equation = e + fM + g 1,,^^ 


considering all earthquakes of magnitude less than 6. 


Province 

e 

Ge 

f 

CTi 

g 


G 

Northern 

1 206 

1.300 

0.731 

0.350 

0.357 

0.270 

0.51 

Northwest 

0.123 

0.605 

0.023 

0.148 

0.963 

0.081 

0.27 

Peninsular 

5.08 

0.878 

-0.603 

0.900 

0.8787 

0.240 

0.69 

Indian sub. 

1 472 

1.000 

0.202 

0.270 

0.704 

0.140 

0.70 







Magnitude ^ Magnitude 







Magnitude 








Magnitude 
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Figure 4.7 Plot of residues of the regression analysis for Northern province considering 
all earthquakes. 
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Figure 4.8 Plot of residues of the regression analysis for Northwest province considering 
all earthquakes. 
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Figure 4.9 Plot of residues of the regression analysis for Northeast province considering 
all earthquakes 
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Figure 4.10 Plot of residues of the regression analysis for Peninsular province 
considering all earthquakes. 






Chapter 5 

Discussion and Conclusions 


5.1 Introduction 

India has experienced a large number of earthquakes in the past. But very few 
strong motion records have been obtained from these earthquakes. On the other hand, the 
intensity isoseismals are available for a substantial number of earthquakes and it is 
pmdent to use these for development of attenuation relationships. 

In this study, 59 earthquakes from different parts of the country have been 
considered. The isoseismals are collected from Memoirs of Geological Suiwey of India 
and some other publications. All the isoseismals along with their intensity area data have 
been compiled. Considering the regional geology and tectonics, and the standard 
deviation obtained from the regression analysis, the Indian subcontinent is divided into 
SIX intensity attenuation province, viz., Northwest, North, Northeast, Peninsular, 
Foredeep, and Subduction Provinces. Apart from these, a single attenuation relationship 
has also been derived for the entire Indian subcontinent treating it as a single region. The 
earthquakes are divided into two categories, one with magnitude greater than or equal to 
6.0 and the other with magnitude less than 6.0. These two categories and their 
combinations have been separately analyzed. 

5.2 Discussion of Results 

5.2.1 Type of Equations and Earthquake Parameters Used 

In this study three different earthquake parameters (M, 4) are used. The ten 

different equations used are described in the Section 4.5.7. Equations 3 (involving 

magnitude and both the linear and logarithmic distance tenn), 6 (involving maximum 

intensity and both the linear and logarithmic distance term), 7 (involving maximum 

intensity, magnitude and both the linear and logarithmic distance tenn) and 10 (involving 
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epicentral intensity and both the linear and logarithmic distance tenii) have two terms, 
cR and d logi„{R), for distance while others have only one of these terms for distance. 

This is to find the predominance of either of the distance terms indicating whether 
logarithmic or linear attenuation with distance dominates a particular region, discussions 
here are made only for Equations 3, 6, 7, and 10 since they involve both distance terms. 

Equation 7 contains both maximum intensity {Imax) and magnitude {M), and fits 
better than equations 3 and 6. Equations 7 and 10 are used to study the regional 
dependency for intensity attenuation. 

5.2,2 Selection of the Best Equation 

Equations 10, involving epicentral intensity (/o), fits best for the available data. 
The estimation of epicentral intensity {Ig) involves a graphical and iterative method 
wherein the value of h, depends on the database of intensity and distance. Hence, I„ is 
really not an independent variable indicating size of earthquake or maximum intensity of 
shaking. On the other hand, magnitude and maximum intensity are independent variables. 
Therefore, it is reasonable to expect that regression analysis will give higher standard 
deviation values when using magnitude or maximum intensity as a variable. 

The estimation of epicentral intensity for future earthquake is not possible since it 
IS a hypothetical intensity. Hence, Equation 10 is not very useful for seismic risk 
estimation studies even though it gives smaller values of standard deviation. Equations 3 
and 6 are better suited for attenuation estimation of future earthquakes as they involve 
magnitude and maximum intensity, which may be assumed for a given seismic risk 
scenario. In the present study, a regression analysis has also been performed to estimate 
epicentral intensity (4) from given magnitude and maximum intensity (Tables 4.17 to 
4.20). The standard deviation in this case is quite high. Hence, Equation 3 or 6 is better 
than Equation 10 for seismic risk estimation studies. 
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5.2.3 Regional Dependency 

The Indian subcontinent is divided into six attenuation provinces. Performing the 

Mest (r-static) on the coefficients of the proposed equations, the possibility of combining 
the different attenuation provinces can be obtained. Table 5.1 shows the possibility of 
combining the different provinces. The t-test performed on the coefficients of Equation 
1 0 suggests that north, northeast and subduction zone earthquakes can be combined into 
one set, and northwest and peninsular earthquakes can be combined into another set. As 
per Equation 7, northwest, north and peninsular earthquakes can be combined into one 
set. Therefore peninsular and northwest region can be combined in either of the cases, but 
Figures 5.1 to 5.3 show that they have different attenuation rate. Moreover large standard 
deviation of northwest earthquakes may have given the positive t-test. Hence in the 
present study these have been treated separately. 

5.2.4 Magnitude Dependency 

Keeping in mind that the smaller earthquakes attenuate more rapidly than the large 
earthquakes, the earthquakes are divided into two groups: those with M<6 and those with 
M><d. The t-test performed on these groups shows that both exhibit different 
characteristics (Table 5.2) and could not be combined. Figures 5.4 to 5.6 show attenuation 
trend of smaller and larger earthquakes for different provinces, both indicating quite 
different attenuation trends. 

5.2.5 Choice of Distance 

In the present study, both hypocentral distance and epicentral distance have been 
used. Hypocentral distance fits better for Equations 3 and 10, while epicentral distance 
fits better for Equation 6. Equation 7 involves both maximum intensity and magnitude 
and hence, hypocentral distance and epicentral distance both fit well for Equation 7 
depending upon the dependence on magnitude and maximum intensity respectively. 
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Equation 6 involves maximum intensity, which includes the effect of focal depth, and 
hence, use of hypocentral distance does not give better results in this case (Table 5.3). 

Plot of intensity versus hypocentral distance, and intensity versus epicentral 
distance, are shown in Figure 5.7. For large distances both curves give similar result. 
Subduction earthquakes show considerable difference as they have large focal depth. 

5.2.6 Effect of Isoseismal Modeling 

The isoseismals can be modeled as a family of concentric circles or ellipses. In 
case of ellipses, relationships for intensity decay have been developed both along the 
major and minor axes of the isoseists. Figure 5.8 shows the average, longitudinal and 
transverse attenuation for different attenuation provinces. There is considerable difference 
in the attenuation rate for longitudinal and transverse directions, as large as one to two 
intensity units. Average attenuation will overestimate the attenuation in the longitudinal 
direction and underestimate that in the transverse direction. Hence, where it is possible to 
estimate the fault direction and therefore, the direction of major axis, one should use the 
longitudinal and transverse attenuation for estimation of seismic risk. 

Subduction zone earthquakes show very large difference in the attenuation rate in 
the longitudinal and transverse directions, as large as two to three intensity units, and 
hence, one should be very careful in using the appropriate equation for the subduction 
region. 

5.2.7 Approximate Decay Pattern 

The attenuation relationships contain two temis for distance, c/? representing 
material absoiption and cl log,g{R) representing geometric spreading. Figures 5.9 to 5.15 
compare predominance of either of these two temis for different provinces. Tables 5.4 to 
5 6 give the standard deviation and R' values for Equations 8, 9 and 10. In all the regions 
except subduction zone, the logarithmic decay term dominates. The subduction zone 
earthquakes are dominated by the linear terra. 
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Figure 5.9 shows the plot of attenuation of intensity with different distance 
terms, c/?, cl log, ^{ r) and both, for different provinces. North, northeast, peninsular, 

foredeep and Indian subcontinent earthquakes show logarithmic decay for first 200, 500, 
450, 600 and 500 tan, respectively, and after that both distance terms dominate equally. 
On the other hand, northwest and subduction earthquakes show more or less logarithmic 
and linear decay, respectively. The results are presented in Table 5.1. 

5.3 Attenuation for Different Provinces 

The Indian subcontinent is divided into six attenuation provinces viz., north. 

Northeast, Northwest, Peninsular, Foredeep and Subduction zone. It is found that 
earthquakes in different province follow different attenuation depending upon the type of 
equation used. Hence attenuation of different provinces is described as per the type of 
equation used. 

5.3.1 Equation Involving Magnitude 

Equation 3 has magnitude as the earthquake parameter. The attenuation of each 
province varies according to the magnitude dependency (constant bi). It is found that the 
earthquakes of northeast province have the maximum magnitude dependency. The value 
of constant bi is close to 2 /.e., an increase in magnitude by one unit can cause an increase 
in intensity by about two units, while the values of constant b} for other provinces are 
around one unit. 

Foredeep and subduction earthquakes exhibit quite different characteristics than 
the other provinces. They show very slow attenuation for first 100 km (Figure 5.16) and 
fast decay after 100 km. The other provinces show more or less similar attenuation rates. 

5.3.2 Equation Involving Maximum Intensity 

Equation 6 has maximum intensity as the earthquake parameter and the 

attenuation of earthquakes in different provinces depends upon the maximum intensity 
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dependency (constant bi). It is found that the subduction and foredeep earthquakes have 
the highest dependency on maximum intensity. 

Subduction zone earthquakes show slow attenuation rate for first 100 km (Figure 
5.17) and then rapid decay rate. Foredeep earthquakes show the slowest decay, while 
northern earthquakes have the fastest decay. 

5.3.3 Equation Involving Epicentral Intensity 

Equation 10 gives drop of intensity from epicentral intensity with distance. It has 

only distance dependency. The attenuation is only due to distance irrespective of the size 
of earthquake. 

Figures 5.1 to 5.3 give the attenuation of different province earthquakes with 

distance. The smaller earthquakes attenuate very fast and earthquakes in northwest 

province have fastest decay while peninsular earthquakes have smallest decay. Among 

larger earthquakes, subduction zone earthquakes have slowest decay for first 400 km and 

after that northeast earthquakes have slowest decay. Northern earthquakes have fastest 

decay for first 200 km and after that peninsular earthquakes have fastest decay. 

5.4 Comparison of Attenuation Pattern of Whole India Subcontinent with United 
States 

For the seismic risk analysis, the attenuation relationships developed in United 
States are often used for Indian projects. Hence, it is interesting to compare the intensity 
attenuation trend of Indian earthquakes with those of United States. Figures 5.18 to 5.20 
compare the Indian attenuation relationships with that of Unites States, and it is found 
that the intensity attenuation pattern of Indian earthquakes matches very well with that of 
San Andreas province of United States. Hence, it is suggested that the strong motion 
attenuation relationships could be used for Indian earthquakes. 
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5.5 Comparison with the Bhui 2001 Earthquake (M 1.9, Imax 10) 

Equations 3, 6, 7 and 1 0 of peninsular province and Indian subcontinent have been 

compared with the recent Bhuj earthquake of January 26, 2001 (Figures 5.21 to 5.24). It 
is found that Equations 10 (involving epicentral intensity and both the linear and 
logarithmic distance tenn) for Indian subcontinent fits best. The Equations 3, 6 and 7 of 
Indian subcontinent fits better than that of peninsular. In any of the case. Equations of 
Indian subcontinent fit better than the peninsular region, the probable cause could be the 
non-availability of intensity data for large earthquakes in peninsular province. 

5.6 Conclusions 

1 . Out of different attenuation Equations, Equation 10 (involving epicentral intensity and 
both the linear and logarithmic distance term) fit best of all the provinces but Equation 
3 (involving magnitude and both the linear and logarithmic distance term) or 6 
(involving maximum intensity and both the linear and logarithmic distance term) can 
be used for seismic risk evaluation. 

2. The northeastern earthquakes have very large dependence on magnitude, and the 
attenuation pattern changes significantly with slight increase in the values of 
magnitude. 

3. Hypocentral distance fits best for Equation 3 (involving magnitude and both the linear 
and logarithmic distance term) and Equation 10 (involving epicentral intensity and 
both the linear and logarithmic distance term), while epicentral distance fits best for 
Equation 6 (involving maximum intensity and both the linear and logarithmic distance 

teim). 

4. Subduction and foredeep earthquakes show different characteristics than the rest. 
They have relatively high epicentral tract than that of others indicating maximum 
damage areas. Subduction earthquakes have linear attenuation, while the other 
provinces follow a logarithmic decay. Northern, northeastern, peninsular and foredeep 
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earthquakes show logarithmic decay for first 200, 500, 450 and 600 km respectively, 
after which both distance term dominates equally. 

5. Small intensity earthquakes attenuate faster than large intensity earthquake, and 
hence, can be separated out. 

6. Among the smaller earthquakes, northwestern earthquakes attenuate fastest, while the 
peninsular earthquakes have the slowest decay rate. 

7. It is found that the attenuation along the transverse direction of the isoseismal is faster 
than that along the longitudinal direction. The difference is as high as one to two 
intensity units. Assuming equivalent circles for isoseismals introduce approximation, 
and hence, separate attenuation relationships have been developed along the minor 
and major axis of the isoseismals. 

8. The attenuation trend of India matches very well with that of San Andreas region and 
hence it is suggested that the strong motion attenuation relationships of San Andreas 
could be used for Indian earthquakes. 

5.7 Scope for Future Studies 

1 . Equations involving magnitude (and maximum intensity) as the earthquake parameter, 
have dual dependency, both on magnitude (and maximum intensity) and distance. 
Hence, it is not possible to get intensity decay independently as a function of distance. 
One can develop relationships separating their dependency. 

2. For the present study the method of least square is used in the regression analysis. The 
two-stage regression method or the method of maximum likelihood, which are 


reported to be most robust, can be used. 
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Table 5.1 Possibility of combining the regions as per t-test values for the comparison of 
constants of equation 10 for hypocentral distance (in bracket are those of 
equation 7). 


Province 

North 

Northeast 

Northwest 

Peninsular 

Fooredeep 

Subduction 

North 






||■|[|| 

Northeast 







North vrest 


X(X) 





Peninsular 

A((/) 

XiX) 

✓(✓) 




Foredeep 

X{X) 

X(X) 

X{X) 

^iX) 



Subduction 


^iX) 

X(X) 

XiX) 

XiX) 



Table 5.2 Possibility of combining the smaller and larger earthquakes group, as per t- 
test values for the comparison of constants, of equation 10 for hypocentral 
distance (in bracket are those of equation 7). 


Province 

North 

Northwest 

Peninsular 

North 




Northwest 

XiX) 



Peninsular 

XiX) 

XiX) 



Table 5.3 Best distance parameter for different type of equations. 



Best distance parameter 

1,2,3(M) 

Hypocentral 

4, 5, 6 (Imax) 

Epicentral 

7 (M, In, ax) 

Hypocentral or Epicentral 

8, 9, 10 do) 

Hypocentral 
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Table 5.6 Table showing standard deviations and values for different provinces as 

per equation 8, 9 and 10 respectively, for transverse attenuation. 


Province 

R(Eq-8) 

Log(R)(Eq-9) 

With Both (Eq-10) 

S.D. 


S.D. 

R^ 

S.D. 

R^ 

North 

1,037 

0.610 

0.769 

0.785 

0.773 

0.789 

Northwest 

1.170 

0.815 

0.758 

0.922 

0.696 

0.935 

Northeast 

1.107 

0.803 

0.495 

0.960 

0.299 

0.986 

Peninsular 

0.976 

0.780 

0.581 

0.922 

0.403 

0.963 

Subduction 

0.674 

0.919 

0.778 

0.892 

0.545 

0.951 

Foredeep 

1.170 

0.794 

0.653 

0.936 

0.551 

0.956 

North (M<6) 

0.82 

0.804 

0.804 

0.811 

0.764 

0.838 

Northwest(M<6) 

0.628 

0.758 

0.597 

0.781 

0.588 

0.796 

Peninsular(M<6) 

0.521 

0.832 

0.441 

0.879 

0.396 

0.908 

North (M>6) 

0.782 

0.687 

0.533 

0.854 

0.524 

0.872 

North west(M>6) 

1.101 

0.867 

0.737 

0.929 

0.613 

0.952 

Peninsular(M>6) 

1.083 

0.790 

0.431 

0.967 

0.370 

0.976 


Table 5 . 7 General trend of attenuation pattern for different attenuation provinces. 


SI. No 

Province 

Average 

Longitudinal 

Transverse 

1 

North 

Logarithmic 

Logarithmic 

Logarithmic 

2 

Northwest 

Logarithmic 

Logarithmic 

Logarithmic 

3 

Northeast 

Logarithmic 

Logarithmic 

Logarithmic 

4 

Peninsular 

Logarithmic 

Logarithmic 

Logarithmic 

5 

Foredeep 

Logarithmic 

Logarithmic 

Logarithmic 

6 

Subduction 

Linear 

Linear 

Linear 

7 

North(M<6) 

Both 

Both 

Both 

8 

Northwest(M<6) 

Logarithmic 

Logarithmic 

Logarithmic 

9 

Peninsular(M<6) 

Logarithmic 

Logarithmic 

Logarithmic 

10 

North(M>6) 

Logarithmic 

Logarithmic 

Logarithmic 

11 

Northwest(M>6) 

Logzirithmic 

Logarithmic 

Logarithmic 

12 

Peninsular(M>6) 

Logarithm 

Logarithmic 

Logarithmic 
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Table 5.8 Approximate pattern of type of attenuation patterns with distance. 


SI. No 

Province 

Smaller 

distance 

Larger 

distance 

1 

North 

< 200km Log 

> 200 Both 

2 

Northwest 

Logarithmic 

Logarithmic 

3 

Northeast 

< 500km Log 

> 500 Both 

4 

Peninsular 

< 450km Log 

> 450 Both 

5 

Foredeep 

< 600km Log 

> 600 Both 

6 

Subduction 

Linear 

Linear 

7 

North(M<6) 

< 120km Log 

> 120 Both 

8 

Northwest(M<6) 

< 90km Both 

> 90 Both 

9 

Peninsular(M<6) 

< 1201an Log 

> 120 Both 

10 

North(M>6) 

Logarithmic 

Logarithmic 

11 

Northwest(M>6) 

Logarithmic 

Logarithmic 

12 

Peninsular(M>6) 

Logarithm 

Logarithmic 
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Figure 5.1 Attenuation of intensity with hypocentral distance for different provinces 
using equation 10 considering all earthquakes. 



Figure 5.2 Attenuation of intensity with hypocentral distance for different provinces 
using equation 10 considering earthquakes of magnitude greater than or 
equal to 6. 
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Figure 5.3 Attenuation of intensity with hypocentral distance for different provinces 
using equation iO considering earthquakes of magnitude less than 6. 
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Figure 5.4 Plot of attenuation rates of different earthquake group for northern province 
as per equation 10. 



Figure 5.5 Plot of attenuation rates of different earthquake group for northwest province 
as per equation 10. 



Figure 5.6 Plot of attenuation rates of different earthquake group for peninsular province 
as per equation 10. 
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Figure 5.8 Plot of attenuation along average, longitudinal and transverse direction for 
different provinces as per equation 10. 
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Figure 5.9 Attenuation pattern for different attenuation provinces with distance as 
parameter (equations 8, 9 and 10) 
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Figure 5.10 Plot of attenuation with different distance parameters for north province 
considering earthquakes of M<6 (equations 8, 9 and 1 0). 



Figure 5.11 Plot of attenuation with different distance parameters for northwest province 
considering earthquakes of M<6 (equations 8, 9 and 10). 



Figure 5.12 Plot of attenuation with different distance parameters for peninsular 
province considering earthquakes of M<6 (equations 8, 9 and 10). 
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Figure 5.13 Plot of attenuation with different distance parameters for north province 
considering earthquakes of M>6 (equations 8, 9 and 10). 



Figure 5.14 Plot of attenuation with different distance parameters for northwest 
province considering earthquakes of M>6 (equations 8, 9 and 10). 



Figure 5.15 Plot of attenuation with different distance parameters for peninsular 
province considering earthquakes of M>6 (equations 8, 9 and 10). 
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Figure 5.16 Plot of average attenuation for different provinces for magnitude 8.5, 7.0 
and 5.5 using equation 3. 
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Figure 5.17 Plot of average attenuation for different provinces for maximum intensity 
11,9 and 6. 
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5.18 Comparison of attenuation pattern of India to that of United States as given 
by Chandra (1979). 



Figure 5.19 Comparison of attenuation pattern of India to that of others as given by 
Gupta and Nuttli (1976). 
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Figure 5.20 Comparison of intensity patterns of India with that of United States 
(Chandra 1982). 
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Figure 5.21 Comparison of Bhuj 2001 earthquake data with that of derived equation for 
Peninsular province (Equations 3, 6, and 7). 



Epicentral Distance (km) 


Figure 5.22 Comparison of Bhuj 2001 earthquake data with that of derived equation for 
Indian Subcontinent (equations 3, 6, and 7). 
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5.23 Comparison of Bhuj 2001 earthquake data with that of derived equation for 
Indian Subcontinent (equation 10). 



Figure 5.24 Comparison of Bhuj 2001 earthquake data with that of derived equation for 
Peninsular region (equation 10). 
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